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This work reviews the basic theoretical aspects, the main observational evidences and the recent apphcations of 
gravitational lensing in the Universe. The article is aimed particularly at providing the readers who don't work on 
gravitational lensing a relatively easy introduction to this active research field in today's astrophysics. 
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INTRODUCTION 

One of the most important tasks for astronomers and physicists is to study the matter distribution in the 
Universe. Based on the assumption of "light-traces-mass" , the map of the Universe can be directly drawn 
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from the measurements of the apparent positions of various luminous objects on the sky (two dimensions) 
and their distances from the earth (one dimension). This method has been widely used nowadays and has 
revealed the existence of large-scale structures, such as the "Great Wall", voids, filaments, etc. However, 
the dark matter puzzle in today's physics and astrophysics casts doubt on the hypothesis of using luminous 
objects as the tracers of the total matter distribution of the Universe. Indeed, astronomical observations can 
only give rise to the distributions of those celestial objects that have electromagnetic radiation strong enough 
to be captured by telescopes, which might not reflect at all the real matter distribution in the Universe. 

Dynamical methods are traditionally used for the determination of the masses gravitationally bound 
in the celestial bodies, which have successfully led to the discoveries of the excess of dynamical masses in 
galaxies and clusters of galaxies as compared with their luminous masses. Nevertheless, the employment of 
dynamical analysis requires that the systems are in the state of dynamical equilibrium, while it has remained 
unclear to date whether the large gravitationally-dominated systems like clusters of galaxies have reached 
the virialized cosmo-dynaniic;al state, especially at high redshifts. Furthermore, the mass determination 
using dynamical method relies on the detailed knowledge of the celestial bodies: The rotational velocities 
and/or velocity dispersions should be well measured in galaxies and galaxy clusters in order to estimate their 
dynamical masses. Yet, this turns out to be quite difficult for the distant galaxies and galaxy clusters. 

A new method of mapping the total matter distribution (luminous+dark) in the Universe stems from 
the effect of gravitational lensing, which has been available only for about 16 years since the discovery of the 
first gravitationally Icnsed double quasar 0957+561 A, B (Walsh, Carswell and Weymann, 1979). The masses 
derived from gravitational lensing reflect the total matter contained in the lensing objects, independent of 
whether or not the lensing systems have reached the virial equilibrium. Therefore, gravitational lensing 
provides an independent way to test the evolution and to set constraints on the possible form of matter 
distributions of the lensing objects. Moreover, studies of gravitational lensing open a possibility of weighing 
the unseen matter in the Universe by comparing the gravitational masses deduced from lensing with the 
luminous masses estimated from optical observations. Furthermore, the effect of gravitational lensing mag- 
nifies the apparent luminosities of background objects, making the intrinsically faint sources enter into the 
detection thresholds of telescopes, which acts in fact like a "gravitational telescope" . 

Recall that modern cosmology is based on the so-called cosmological principled , which assumes a spatially 
isotropic and homogeneous matter distribution of the Universe. The isotropy of the Universe has been well 
demonstrated by the measurement of the microwave background radiation (Smoot et al., 1992), whilst 
the homogeneity turns to be somewhat hard to describe quantitatively. Actually, it is not so clear on what 
scales the Universe can be treated as homogeneous, though the largest coherent structure seems to have 
scale of ~ 100 Mpc. Matter condensations occur on scales up to a few tens of megaparsecs: planets, stars, 
galaxies, galaxy clusters, superclusters, voids. Great Walls, etc. All these matter clumps may affect the 
propagation of light through the effect of gravitational lensing, according to the prediction of the theory of 
general relativity. Therefore, gravitational lensing is a common phenomenon in astronomical observations 
and reasonable caution should be exercised in the identification of various celestial bodies. For instance, 
some close double or multiple images may be due to single sources (e.g. quasars) gravitationally split by the 
intervening objects (e.g. galaxies), and the arclike images could be the result of the gravitationally distorted 
background galaxies by the foreground galaxy clusters. 

The history of studies of gravitational lensing can be divided into three periods: 1704 - 1964, 1964 - 1979 
and 1979 - present. Newton (1704) addressed the question nearly three hundreds years ago if the celestial 
bodies could bend light rays. The deflection of light by a spherical body of mass M based on the Newtonian 
mechanics was computed to be a = 2GM/c^^, assuming an impact distance of ^, i.e., the shortest distance 
from M to the light path, c and G are the speed of light and the gravitational constant, respectively. This 
formula was found to imderestimate the deflection angle by a factor of 2 by Einstein in 1915 in terms of his 
gravitational theory. In particular, utilizing the result of general relativity to the Sun predicts a deflection 
angle of 1.75" for a light ray passing near the solar limb. This prediction was very soon confirmed during 
the solar eclipse in 1919 by a team led by Eddington (Dyson, Eddington and Davidson, 1920). Although 
some progress had been made since then on the theoretical aspects of light bending, not much interest had 
been really drawn on this field until 1963-1964 when the geometry of lensing was studied independently by 
Klimov (1963), Liebes (1964) and Refsdal (1964a). Their work can be considered pioneering in the sense that 
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they set up the foundation of the modern theoretical research on gravitational lensing. During the period of 
1964 - 1979 some important progress was made in the computations of deflection angle (e.g., Bourassa and 
Kantowski, 1975; 1976) and the light propagation in the model of an inhomogeneous Universe (Press and 
Gunn, 1973; Dyer and Roeder, 1972; 1973; 1974). These achievements have played an important role in the 
lensing studies after the detection of the first multiple images of quasar 0957+561A,B by Walsh, Carswell 
and Weymann in 1979. In the decade following this landmark, the detections of new lensing phenomena have 
increased dramatically, including multiply- imaged quasars, giant luminous arcs and arclets, radio Einstein 
rings, microlensing events, whilst the theoretical investigation of gravitational lensing has concerned many 
interesting subjects of modern cosmic physics, such as determinations of total masses of the lensing systems, 
determinations of Hq and f^o, explanation of the associations of bac;kgroun(l sources with foreground objects, 
searches for dark matter candidates, etc. Today, gravitational lensing has become one of the most active 
fields in cosmological research. 

It is surely impossible to cover all the topics of lensing research in this review as it has expanded very 
rapidly in the past years. The purpose of this work is to concentrate on the basic theories and the new 
progress of gravitational lensing as well as their applications in astrophysics. Alternatively, this review will 
not trace the history of development of gravitational lensing but follow the scales of lenses from compact 
objects to large-scale inhomogeneities. 



1 BASIC THEORY 
1.1 Deflection of light 

Gravitational lensing is based on the theory of general relativity, which predicts that light rays would be 
bent when they pass near a massive body. For a pointlike mass M the deflection angle of light rays is 

«=^. (1) 

where ^ represents the impact distance of light rays. This deflection angle is twice as large as the value 
predicted from Newtonian mechanics. 

For an extended mass distribution, the deflection of light passing through the mass system cannot be 
simply obtained in the frame of general relativity. Only for some special matter distributions can the metric 
and the solution to the photon geodesic equations, hence the angle of light bending, be found exactly. For 
instance, in a spherical uniform matter distribution the deflection of light, to first order in the Newtonian 
gravitational potential, is (Wu, 1989a) 



4GM 



^2 \ 3/2- 



(2) 



Here M and R are the total mass and the radius of the massive sphere, respectively. Another matter profile 
that has been commonly used for galaxies and clusters of galaxies is the singular isothermal sphere model 
having mass density 

where cr„ measures the line-of-sight velocity dispersion. The deflection angle of light rays can be found to 
first order in {ayjcf to be (Wu, 1989a) 

a = 47r4, (4) 
i.e., a constant deflection independent of the impact distance. 
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For photons traveling in the gravitational field of an irregular matter distribution, the deflection angle of 
light rays is often computed in a linearized Einstein approximation 



a 



2 f°° 

- / V<l>dt, (5) 

^ J — oo 



in which (j) represents the Newtonian gravitational potential of the matter distribution. An equivalent but 
simple formula is obtained by dividing the system into a number of pointlike masses (Mj) and then summing 
up the contributions from each small mass piece 

i ' ' 

where r — ri is the impact vector of light rays from the mass unit Mj. The integral form is 

« = f / jl^^n^d^'- (7) 

Here is the surface mass density at position ^' obtained by projecting all the mass along the line of 

sight onto the "lens plane". In particular, for a spherical matter distribution p(r) the above equation can 
be simplified to be 

4Gm(g) 

where m(^) is the total projected mass along the light-of-sight enclosed within the impact distance ^. If p{r) 
is confined within the radius R, then 

m(C) =M-4n y/^i^^^rp{r)dr, ^ < R; .g. 

m(C) = M, ^>R, 

in which M is the total mass of the spherical system. Therefore, if the deflector has spherical symmetry, 
the light bending can be obtained simply by replacing the pointlike mass M of eq.(l) by the projected mass 
m(^). For example, replacing the mass density in eq.(9) by a constant and inserting m(^) for ^ < R into 
eq.(8) recover the result of eq.(2). It should be pointed out that eq.(9) holds true even if the deflector has 
no boundary or infinite mass, in which 

/oo 
(r - - ^2)rp{r)dr. (10) 

So, the deflection angle [eq.(4)] of a singular isothermal sphere can be easily obtained by inserting the density 
proflle eq.(3) into eqs.(8) and (10). 



1.2 Lensing geometry and Vcirious lenses 

Suppose that a light-ray from a distant source at redshift Zs passes through or near a massive system at 
redshift Zd (zd < Zg) and then reaches the observer at redshift z = (Figure 1). If there wore no intervening 
massive system, the light from the source would have arrived at the observer along a straight line. The 
gravitational field of the massive system now bends the light ray, causing the direction of light to be changed 
by an angle of a. If we use 77 (/3 in angle) to denote the true position of the source or the alignment 
parameter in the source plane and ^ {G in angle), the observed positions of the images in the lens plane, we 
have the following geometrical relation, namely, the "lensing equation" : 

^4 - 77 = Ddsoc, (11) 
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Figure 1: Scheme of basic geometry of gravitational Icnsing. The true position of the distant source is /? 
with an angular diameter distance Dg from the observer. The light rays with an impact distance 9 would 
be deflected by an angle a by the gravitational field of an intervening massive object, resulting in multiple 
and distorted images. 



or 



6-/3 



D. 



-a 



(12) 



where D^, Dg and D^g are the angular diameter distances to the deflector, to the source and from the 
deflector to the source, respectively. As seen from eq.(7), the deflection angle a can be written as the 
gradient of a two-dimensional potential ip 

= (13) 



and 



in which the quantity Sc is called the critical surface mass density 

Dg 



(14) 



(15) 



and the physical surface mass density term S should satisfy the two dimensional Poisson's equation (Bland- 
ford and Narayan, 1986) 

,2 _ _ (16) 



For a given source position l3, the lensing equation may have several solutions for tlic' image position 
6. As a consequence, we can observe multiple images of a single source on the sky. Furthermore, the 
lensing equation (12) describes the distortion of the surface brightness of background sources according to 
the mapping from source plane to lens plane: 



7'(0)=7(/3)=7(0-V^), 



(17) 



where /' and / are, respectively, the observed and intrinsic surface brightness patterns. The magnification 
factor n describes the change of apparent luminosity of the source, which is characterized by the Jacobian 
for the mapping /3 — > 6 (Schneider, Ehlers and Falco, 1992): 





-1 







here 



K = 

7 = 

71 = 

72 = 



^{ipu +1P22); 



(18) 



V7i +71; 

5(V'll - V'22); 

V'12 = "021, 



(19) 



in which tpij = d^tp/d9id6j. When the line of sight completely misses the defiector, S(0) = and the k term 
vanishes in terms of eq.(16). So, k represents the amplitude of the convergence due to the matter within 
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Figure 2: Image configurations by a black hole (Figure courtesy of C. Zahn and H. Ruder) 



the light-ray (also referred to as Ricci focusing), while the 7 term is the amplitude of the shear due to the 
matter outside the beam (also referred to as Weyl focusing). The latter can be easily calculated in the case 
of spherical deflector, which reads 

AG[m{e)-mm Das 

^ = — — 

and rn{9) = -kO'^ D'^Yj{9) is the mass of the cylinder of radius 6 with a uniform surface mass density equal 
to S(6'). When a beam of light-rays pass through the center of deflector, m(0) = m(0) = and the 7 term 
then becomes zero. 

(l)Pointlike mass as deflector. This model can be considered to be a good approximation for many celes- 
tial bodies like "Jupiters" , stars, black holes, and even galaxies, when the light rays from background sources 
pass outside the deflectors. Solving the lensing equation, using the deflection of eq.(l), yields 



where 



^± = ^(/3±V/3' + 44), (21) 
fAGM Dae 



(22) 



DdDs 

is the critical radius (or aE — OsDa in linear size) corresponding to a ring-like image of a background source 
when (3 — Q, which is often called the "Einstein ring", named after the pioneering work by Einstein (1936). 
6*+ and 9- describe the positions of the two images produced by a point mass M. Their magnifications are 
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1 - {9e/9±\ 



(23) 



If the separation between the two images is too small to be resolved by modern telescopes, as for the 
microlensing events (section 2), the total magnification is often used for their combined effect 

•"^ + 2 

II = 11+ + 11- = ^== (24) 

Here u is defined as u = /3/9e- u = 1 ot (3 = 6e is often taken to be a typical case that characterizes the 
efficiency of the lens, which corresponds to (/i+,/i_) ~ (1.17,0.17) and /i = 1.34 or Am = 0.32 in apparent 
magnitude. Pointlike masses as lenses play an important role in the study of microlensing (section 2.1) and 
in searches for dark matter candidates in the Galactic halo using microlensing effect (section 2.2). 

The image configurations by a pointlike mass is well illustrated in Figure 2 by changing the relative 
positions between the source which is chosen to be the Einstein portrait and the lens which is assumed 
to be a massive black hole. When the source approaches the lens, the two images are elongated, showing 
two arclike structures. Finally, two arc images merge into an "Einstein ring" when the lens is completely 
coincident with the source. 



(2)A singular isothermal sphere. The studies of the flat rotation curves of galaxies and the galaxy/gas 
distributions in clusters of galaxies suggest that the total matter profiles in these systems follow very well 
the singular isothermal sphere (SIS) model [see eq.(3)]. The constant defiection of a = 47r(cr^/c)^ in the 
lensing equation gives 

9±=9e±P (25) 
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Figure 3: Solution to Icnsing equation by ISC. The solid curves represent the lensing equation /3o = 
9q — D[{1 + O^Y^"^ — 1\/9q. The intersecting points with the line (3q =constant give the number and positions 
of the lensed images. 



for P < Oe, where the critical radius is 



= 47r 



(26) 



In particular, the image separation is just the diameter of the Einstein ring: = 6+ + 9- = 26e- Note 
that if the alignment parameter /3 is larger than the Einstein ring, only one image appears instead of two. 
The magnifications of the images are simply 



It 



9e 
0± 



l±47r- 



(27) 



(3)A softened singular isothermal sphere. SIS, though simple, is an unphysical model because the mass 
density reaches infinity at the center. Instead, a SIS with a finite core radius of Tc (or Oc in angle) (ISC) 
seems to be more reasonable for the matter distributions of galaxies and galaxy clusters (Hinshaw and 
Krauss, 1987): 



p{r) = 



1 



(28) 



This density profile reduces to SIS when rc 
mass within ^ are 

m 



2-kG r'^ + rl' 

or r Tc- The surface mass density and the total projected 



(29) 



m(e) = ^(v/F+^- 

respectively. The deflection of light can be directly obtained from eq.(9) and then, the lensing equation reads 

1 



eo = po + D 



^0 



(30) 



in which Oq and (3q are in unit of Oc- The lensing parameter, defined a& D = {AiraH c^){DdDds/rcDs), 
determines the number of the solutions (Figure 3) (Wu, 1989b). In the case of £> < 2, ISC always produces 
a single image, while it may result in three images for > 2 if /3o is sufficiently small. The intersections 
of the lines /3o =constant with the curves (3o = 0o — D{y/l + 9q — l)/9o give the solutions to the lensing 
equation. The magnification is found to be 



1-D 



01 



l + D 



01 



1 



D 



1 



(31) 



(4)0ther spherical models. Two other models which are also frequently adopted for the matter distri- 
butions in galaxies and in galaxy clusters are the King model (or the modified Hubble model) and the de 
Vaucouleur model (or the r^/^ law). Their lensing properties are very similar to those of ISC and are sum- 
marized in Table 1. 



(5)Asymmetric lenses. In principle, the lensing geometry can be established for any kind of geometrically- 
thin matter inhomogeneity in the Universe through eqs.(7) and (12). Among them, the properties of elliptical 
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Table 1: Gravitational Lensing Models 



Figure 4: Image configurations by an elliptical lens. The source planes are on the left and the corresponding 
images arc on the right. The solid lines are the caustics and the dashed lines are the corresponding critical 
lines, (from Blandford and Narayan, 1992) 



lenses have been thoroughly studied (e.g., Bourassa and Kantowski, 1975; 1976; Blandford and Kochanek, 
1987; Narasimha, Subramanian and Chitre, 1987; Schramm, 1990; Wallington and Narayan, 1993; Kassi- 
ola and Kovner, 1993, etc.). A simple elliptical lens assumes the following "non-singular pseudo-elliptical 
isothermal potential" 



Att- 



c Ds 

where e is the ellipticity and 6c is the core radius of the lens. Solving the lens equation 

13 = e - vi^, 

one has (Blandford and Kochanek, 1987; Wallington and Narayan, 1993) 



Pi 

f32 



gj, (1^ 



(32) 



(33) 



(34) 



where 6e is the Einstein radius in SIS [eq.(26)]. Magnifications of the images can be calculated using eq.(18) 



dOi 802 



862 dOi 



(35) 



Inverting eq.(34) gives the image positions in the lens plane for pointlike sources. An extended source can 
be considered as a set of pointlike sources and its image configuration produced by an elliptical lens can be 
drawn through the above equation for each source element. In principle, the images with arbitrary shapes 
can be produced by an elliptical lens as long as the position and the shape of the extended source are prop- 
erly chosen. Some image configurations produced in this procedure for a circular source are shown in Figure 4. 



(6)Cosmic strings as lenses. Cosmic strings are topological defects that are believed to be created in the 
very early Universe and could act as the seeds of formations of galaxies, galaxy clusters and even large-scale 
structures (Vilenkin, 1981), although there has been no observational evidence to date for the existence of 
such cosmic strings in the Universe. Cosmic strings, if real, may essentially exist in two forms: straight and 
loop strings. Thus, the long-lived cosmic strings are able to cause gravitational lensing effect on background 
sources. 

The exterior gravitational field of a straight string is {G/c^ = 1) (Gott, 1985) 

ds'^ = -dt^ + dr'^ -I- (1 - AjjLsfr^dct)^ -{- dz^. (36) 

Thus, a cosmic string is completely described by its linear mass density jis (~ 10~^). Consider that a light 
ray propagates in the plane perpendicular to the string (Figure 5), one can easily solve the null geodesic 
equation and find the light deflection in <f) direction to be 
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Figure 5: Geometry of a straight cosmic string 



Figure 6: Loop string and arclike image 

aw47r/is = 2.6" (/is/10"*'). (37) 

The Ught rays that propagate along r and z directions are unaffected by the gravitational field of a straight 
string, according to the metric of cq.(36). Two images of a background source with the separation of 87r/is 
but with the same apparent luminosity appear around the string. Note that eq.(36) adopts a cylindrical 
coordinate rather than the spherical one as in SIS although both lenses produce the impact parameter- 
independent deflections. So, the double images of a background source by a straight string can be identified 
relatively easily from their equal luminosities. 

The metric of a long-lived loop string is unknown today. Its lensing properties can be studied only in 
the linearized gravitational approximation. For the simple case of a "face-on" loop string with radius of a 
in the sky, the light from the source behind the string is bent by (Wu, 1989c) 

a=|0' ^<'^' (38) 

i.e., the light rays which pass through the loop remain unaffected while the light rays outside the loop behave 
in the same way as those by a point mass of M = 27r/tsO at the center. Therefore, if the alignment parameter 
77 of a background source is smaller than the loop radius a, we would always expect to see the original source 
through the loop. In particular, when r\ satisfies a — \Xs{p dsl ^ s) < rj < a, the arclike image of the 
lensed source appears outside the loop. Figure 6 shows such an example, in which the source and loop string 
arc 10 kpc and 100 kpc in radii and located at Zg = 1 and Zd = 0.25, respectively. The equivalent mass of 
the string is ~ lO^'^M©. Therefore, a single giant arclike image, instead of multiple arclike images in the 
case of massive spherical deflectors, can be produced by the loop string. 



1.3 Lensing efficiency 

In the above discussion, it appears that lensing magnification of the apparent luminosities of background 
sources may tend towards infinity in some cases. For example, the impact parameters that satisfy the 
condition (1 — k)^ — 7^ = in eq.(18) correspond to images with an infinitely large magnification, whilst the 
same situation occurs for a pointlike lens when the alignment parameter is 77 = 0. These apparent unphysical 
results arise from the hypothesis that the background source is pointlike. 

In general, the total magnification of an extended source with surface bightness 1(771, 772) can be obtained 
by summing up the contribution ii{rii,ri2) of each source element 1(771,772)^771^772 (Bontz, 1979): 

For simplicity, we consider a circular disk source with uniform surface brightness and radius Rs instead of 
the point source approximation, and we further assume a spherical matter distribution for the defiector. In 
this situation, the maximum magnification iimax corresponds to the case where the centers of the source 
and of the lens and the observer lie perfectly on a straight line, i.e, the Einstein ring shows up. So, eq.(39) 
reduces to the ratio of the area of the Einstein ring to the luminous area of the original source. This ratio or 
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(39) 



Figure 7: Maximum magnification and source radius. The background source is assumed to be a luminous 
circular disk with radius Rs at Zs = 2 while the lens is modeled by a point mass for star [eq.(40)] and SIS 
for galaxy and cluster of galaxy [eq.(41)] at Zd = 0.5. 



maximum magnification is not infinite any longer. Actually, fimax measures the efficiencies of gravitational 
lensing by various lenses for the same source. The larger the Umax is, the stronger the lens would be. For a 
pointlike lens M, 

M„^a. = yi + ^^^; (40) 

and for SIS, 

where Og = Rs/Dg is the angular radius of the source. 

To show how efficiently the various lenses act on the background source, the maximum magnification of 
a background circular disk source is illustrated in Figure 7 for three typical lenses in the Universe; stars, 
galaxies and clusters of galaxies (Wu, 1992a). The star is modeled by a point mass with a solar mass Mq, 
while SIS is adopted for galaxies and galaxy clusters whose velocity dispersions are taken to be 200 km/s and 
1500 km/s, respectively. It is concluded from Figure 7: (l)Compact objects like stars as lenses are capable 
of producing significant lensing effect on sources with sizes smaller than ~ 0.01 pc. Therefore, they can 
affect AGNs (~ 10~^ pc), quasars and normal stars; (2)Lenses on scale of galaxies can magnify any sources 
with size smaller than galaxies themselves; (3)Finally, clusters of galaxies are efficient lenses for nearly all 
luminous objects including galaxies of sizes of ~ 10 kpc. 



2 COMPACT OBJECTS AND MICRO-LENSING 

2.1 Microlensing 

For a pointlike lens M and a background source, both at their typical cosmological distances of Zd = 0.5 and 
Zs = I, the angle subtended by the Einstein radius is [eq.(22)] 

Oe « 1.4 xlO"*^ VMr) arcseconds. (42) 

Thus, a compact star-like lens at cosmological distance gives rise to an image splitting of a background 
source of the order of "microarcseconds" . Hence, the terms "microlens" and "microlensing" are used for the 
compact object and its lensing phenomena so as to distinguish it from the "macrolens" like a galaxy which 
results in the image separation of the order of ~ 1 arcsecond according to eq.(42). Nevertheless, it appears 
to be hopeless to resolve the images microlensed by star-like compact objects in the Universe even with the 
present advanced telescopes. 

Because of the too small separation between the multiply microlensed images, the combined magnification 
of the images turns to be the unique feature that arises from microlensing. Unfortunately, this feature cannot 
be straightforwardly used to detect microlensing events since it is impossible to separate the magnification 
effect from the intrinsic luminosity of the source. However, if the source and/or the lens have a transverse 
motion with respect to the line of sight, the lensing magnification would vary with time, leading to the 
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Figure 8: Microlensing event QSO 2237+0305 (Figure courtesy of R. Stabell) 



Figure 9: Light curves of the four components of QSO 2237+0305. The data are taken from Gumming and 
De Robertis (1995) and the error associated with each datum point is typically 0.05 magnitude. 

variation of apparent luminosity of the source. This effect was firstly predicted by Chang and Rcfsdal in 
1979. Since then, a great number of papers have appeared, trying to apply this property for the explanation 
of QSO variabilities and the searches of dark matter candidates of the Galactic halo (e.g.. Young, 1981; 
Ganizares, 1982; Ostriker and Vietri, 1985; Nottale, 1986; Paczyhski, 1986; Schramm ct al., 1994; etc.). 
The timescalc of luminosity variability of a source by a microlens M can be estimated using the time of the 
source crossing the Einstein radius aE with a relative velocity v. T = 2aE/v. (l)For a local Galactic lens, 
e.g. an object in the Galactic halo (Dd ^ 10 kpc) acting as lens and a star of the Large Magellanic Gloud 
(LMG) {Dg — 50 kpc) being the target source: 



where (also hereafter) Hq = 50/i5o km s^^ Mpc^^ is the Hiibble constant. The timescale of luminosity 
variability of a source due to microlensing depends on the mass of the microlens as ~ {M/Mq)^^^. It turns 
out that the masses of the lensing objects can be determined by monitoring the variations of the apparent 
luminosities of some sources if the distances of lenses and sources as well as the transvcirse speed of the lens 
V are known. Note, however, that this procedure is actually restricted by the sensitivities of the microlensing 
observations, which cannot include the events with timescales lasting both shorter than the observing periods 
of sampling and longer than the observing coverage. The present observations are then sensitive to those 
events with timescales ranging from ~ 10 minutes to ~ 10 years. 

The convincing evidence that a distant source is microlensed has been found so far only in QSO 2237+0305 
{zs — 1.695) associated with a foreground spiral galaxy (zd = 0.0394). This lens system was discovered by 
Huchra et al. in 1985, in which the quasar appears to be coincident with the nucleus of the galaxy. The 
subsequent observations (De Robertis and Yee, 1988; Yee, 1988; Schneider et al., 1988) show that this 
quasar actually consists of four components with the maximum separation of 1.8" (see Figure 8), which is 
also referred to as the Einstein cross. The photometric monitoring of QSO 2237+0305 has been made 
for several years and the brightness variations of each component are shown in Figure 9. Overall, the light 
curves of four components show no apparent correlations, indicating that their variations must be due to 
the microlensing magnification by the compact objects and their relative motions in the spiral galaxy. Note 
that the time delay between the components from the microlensing is estimated to be about one day. To be 
specific, a sharp variation of at least 0.2 magnitude in the component A was detected within 26 days during 
1988-1989 (Irwin et al., 1989; Gorrigan et al., 1991), while the other components didn't exhibit a similar 
feature. Another strong microlensing occurred in component B in 1991: the brightness of B is magnified by 
a factor of about 0.5 magnitude (Yee and De Robertis, 1992; Racine, 1992). As a comparison, components 
G and D appear to be relatively stable. 




(43) 



and (2)for a lens at cosmological distance {za = 0.5 and Zs = 1) 




(44) 
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Figure 10: Light curve of a pointlike source passing behind a microlens (marked by "•") with a relative 
velocity of v. The impact distance is 0.25aE and as is shown by the dotted lines. 



2.2 MACHO searches 



The inferred masses from the rotation curves of galaxies, including our Galaxy, are an order of magnitude 
larger than their luminous masses, implying that galaxies are embedded in invisible massive halos (see Ash- 
man 1992 for a recent review). The nature of the dark matter in the halos of galaxies is still unknown today. 
Basically, two kinds of candidates have been proposed: WIMPs (weakly interacting massive particles such 
as axions and supersymmetric neutralinos) and MACHOs (massive astrophysical compact halo objects such 
as brown dwarfs, low mass stars and black holes). WIMPs are some kinds of unknown non-baryonic matter 
which may dominate the Universe according to the standard inflation cosmological model and the Big Bang 
primordial nucleosynthesis, especially if the new measurement of the deuterium abundance in a high rcdshift 
primordial hydrogen cloud (Songalia et al., 1994; Carswell et al., 1994) is confirmed, whilst astronomers 
might favour MACHOs in the sense that MACHOs may cause the observational effect - microlensing. 

If the relative velocity of the microlens M (or the source) transverse to the lino of sight is v, the variability 
of the total magnification of a pointlike source follows eq.(24). Figure 10 shows the light curve of a back- 
ground source passing near a microlens with impact parameter of 0.25aE- The main characteristic features 
of the microlensing light curve arc the achromaticity and timc-symmctry around the point of maximum 
magnification, which can then be distinguished from the known variable star phenomena. 

The probability that a source is gravitationally lensed is described by the so-called optical depth (r). The 
optical depth to microlensing without involving the lenses at cosmological distance is simply the number of 
microlenses or MACHOs inside the "microlensing tube" which has a cross-section of ttq^: 



where n{D(i) is the number density of the MACHOs at distance Dd- Note that within the microlensing tube, 
/U > 1.34 or I Am I > 0.32. Suppose that the Galactic halo is composed of MACHOs, their density profile 
can be estimated from the rotational velocity vq which is about 220 km/s at the position of the Sun. The 
original work by Paczynski (1986) used a SIS for the massive halo, and the subsequent work by Griest (1991) 
modified SIS by introducing a core radius (ISC) which actually does not provide any significant difference 
from SIS in the calculation of r. Employing ISC for the MACHO distribution of the Galactic halo in eq.(45) 
yields 



where r,, is the core radius of ISC, a is the angle between the line of sight to the source and the direction to 
the Galactic center, and Dh is the extent of the Galactic halo. All the distances are measured in unit of RgCi 
the distance to the Galactic center [k, 8.5 kpc). so that x = Dd/RaCi Xs = Ds/Rgc and Xc = Tc/Dgc- 

Now consider a star in our neighbour galaxy, the LMC, as the target. The light rays from a star of 
LMC reach the observer by passing through the LMC halo, the LMC disk and finally, the Galactic halo. 
MACHOs of both the LMC halo / disk and the Galactic halo are able to gravitationally magnify the apparent 
luminosity of the star in the LMC. The contribution of the LMC halo to the microlensing optical depth can 
be computed in a similar way to that of the Galactic halo 




(45) 




(46) 



and 




1 - ^ arctan ^ ' 




(47) 
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Figure 11: Optical depth to niierolensing of the stars of the LMC by the Galactic halo, the LMC halo and the 
LMC disk. The extent of the LMC halo is taken to be 15° (Schommer et al, 1992) and the LMC rotational 
velocity and Vm are both taken to be 79 km/s. 



and 

To = 



^VLMC_Y - 

V c y 1 - arctan ' 

Rh To 



Here Vc is the core radius of the LMC halo described by ISC, /3 is the angle between the line of sight to the star 
and the center of the LMC, and Rh is the extent of the LMC halo. The distances in eq.(47) are measured in 
unit of Dc, the distance to the LMC center (« 50.6 kpc), so that x = D^/Dc, Xg = Dg/Dc and Xc = Tc/Dc- 
The LMC disk can be modeled by an isothermal self-gravitating disk having density (van der Kruit and 
Searle, 1981) n = no exp(— i?//i)sech^(z/zo)7 with the scale length h in the radial direction and the scale 
height zq in the z-direction. The central number density no relates with the maximum rotational velocity 
{vm) through (Freeman, 1970) no = (l/47rMG/izo)(wm/0.62)2, in which M is the mass of the MACHO. 
Furthermore, the inclination of the LMC disk can be simply taken to be 0°, i.e., a face-on distribution of 
the disk matter. The optical depth to microlensing by the disk is then 

t^'xiXs—x) xmnH ,o/l — xcos/3\ , 

TLMCdisk = TO / '-e W^sech' . , / dx, 48 

Jo Xs V {zo/Dc) ) 

and 

° \0.62cJ \hzo)' 

The contributions of the Galactic halo, the LMC halo and the LMC disk to the microlensing optical 
depth for the stars of the LMC are plotted in Figure 11. The halo of our Galaxy gives rise to a nearly 
constant optical depth crossing the LMC disk, tg = 5 x 10"'', while the LMC halo/disk provide an optical 
depth depending sharply on the positions of the stars, which arises from the fact that only the foregroimd 
MACHOs of the LMC are able to act as lenses for the stars of the LMC itself. For the stars near the LMC 
center (/3 < 0°.5), rLMC,haio > 2 x 10"'' and rLMC,disk > 3 x 10"^. Therefore, if the halo of our Galaxy and 
of the LMC are composed of MACHOs, several million stars should be monitored for the discovery of the 
microlensing events in the LMC. 

Following the proposal of Paczyriski (1986), the EROS (Experience de Recherche d'Objets Sombres) 
and the MACHO collaboration commenced in 1990 their searches for microlensing events of the LMC by 
monitoring the brightness of a few million stars in the LMC, and the OGLE (Optical Gravitational lensing 
Experiment) began in 1992 to conduct a similar search in the direction of the Galactic bulge. Three groups 
announced their discoveries at almost the same time in the autumn of 1993: Three events were detected 
in the LMC (Alcock et al. 1993; Aubourg et al. 1993) and six were seen in the Galactic bulge (Udalski et 
al., 1993). During the writing process of this article, the total "local" microlensing events have grown to 
~ 70. These include more than 10 events found by the OGLE collaboration (Udalski et al., 1994a; 1994b; 
1995) and more than 40 events by the MACHO collaboration (Bennett et al., 1994; Alcock et al., 1995a) in 
the Galactic bulge, and 3 events by the MACHO collaboration (Alcock et al., 1995b) and 2 events by the 
EROS team in the LMC [Note that the EROS n°2 candidate may be an eclipsing binary system rather than 
a microlensing event (Ansari et al., 1995)]. Figure 12 shows the light curves of, and the lensing model fits to 
the MACHO microlensing events and Table 2 summarizes the properties of the 5 microlensing candidates of 
the LMC, which can be regarded as the representatives of all the reported microlensing candidates. 

If these events are indeed generated by the microlensing of the MACHOs along the light of sight rather 
than a new kind of variable stars, one can estimate the mass of the MACHOs using the event duration T (i.e., 
the Einstein ring crossing time) and the maximum magnification Umax- At the point where the background 
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Figure 12: The observed light curves of the three MACHO rtiicrolensing candidates and the best-fit theoretical 
lensing model. [From Alcock et al. (1995b)] 



Table 2: Properties of the microlensing candidates of the LMC 



candidate 


RA(2000) 


Dcc(2000) 


media magnitude(V) 


f-^max 


T (days) 


MACHO n°l 


05 15 44.5 


-68 48 00 


19.6 


7.20 


34.8 


MACHO n°2 


05 22 57.0 


-70 33 14 


20.7 


1.99 


19.8 


MACHO n°3 


05 29 37.4 


-70 06 01 


19.4 


1.52 


28.2 


EROS n°l 


05 26 36 


-70 5737 


19.0 


2.5 


27 


EROS n°2 


05 06 06 


-65 58 34 


19.3 


3.0 


30 



star enters into the microlensing tube (denoted by subscript ''min"), eq.(24) reads 



l^min = """"T^ ^ . = 1-34, (49) 



u 



mm ^/u^nin + ^ 



in which Umin = Ij while at the maximum magnification (denoted by subscript "t") 



1^^ = — 7^=7- (50) 



Alternatively, ut is related with Umin through geometrical relation 

/ vT \ 



2a 



where v is the relative velocity of the star or the microlens. These three equations give rise to the mass of 
microlens 

If the MACHOs of Galactic halo are responsible for the observed events, we can take Dd = 10 kpc and 
V = 220 km/s for a numerical estimate. Utilizing the observed microlensing event duration of typically 
T ~ 30 days and the maximum magnification of a few leads to M ^ O.IMq, i.e., sub-solar objects in the 
Galactic halo are likely to be the deflectors for the microlensing events of the LMC. 

However, one cannot conclude from the presently detected microlensing events in the LMC that the 
halos of the galaxies are dominated by ^ 0.1 solar mass objects. In fact, the positions and relative 
velocities v of traverse motion of the lensing objects are two unknown factors in the determination of masses 
of MACHOs. Unless a statistical sample of microlensing events in the LMC is completed, it is in principle 
impossible to draw a decisive conclusion about the masses of the MACHOs in the Galactic halo. Gould 
(1994) argued that the observed optical depth toward the LMC center in the MACHO collaboration is only 
7 — 9 X 10^®, much less than that expected from the Galactic halo made of MACHOs. Using the fact that 
3 events were detected among 9.5 million monitored stars in LMC for 1.1 years, the MACHO Collaboration 
(Alcock et al., 1995a, c) has recently reached a similar microlensing optical depth of 8.8^g x 10~*, nearly 
an order of magnitude lower than the expected optical depth of ~ 5 x 10^^. So, the halo of our Galaxy 
may have not been detected at all. If so, the microlensing events seen by EROS and MACHO collaboration 
may have arisen from the stars of the LMC disk (Wu, 1994b; Sahu, 1994) and self-lensing by a stellar disk 
remains to be an interesting model for further investigation (Gould, 1995). 
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Figure 13: (a) Magnification patterns on the source plane by an ensemble of N (left) and N — 1 (right) 
star-like microlenses and (6)four horizontal (left) and four vertical (right) light curves for a source crossing 
the tracks marked with think black lines in (a). The thick and thin lines correspond to (a)-left and (6)-right, 
respectively. (Figure courtesy of J. Wambsganss) 

2.3 Simulations and observations 

Cosmological compact objects either bounded in galaxies or distributed randomly in the Universe arc capable 
of magnifying temporarily the background sources like quasars, AGNs, etc., resulting in variations of their 
apparent luminosities. Besides the significant difference of timescales between the cosmological microlenses 
and the local ones like those in the Galactic halo and in the LMC [see eqs.(43) and (44)] which we have 
discussed in the above subsection, the optical depth to microlensing arising from the lenses at cosmological 
distance may be a few orders of magnitude larger than the local optical depth, depending on the content 
of compact objects of the Universe. For example, the optical depth to microlensing for a distant source at 
redshift = 3 ~ 4 can be of order of unity if the Universe is composed of compact objects (see Figure 
17). So, one now needs to deal with the problem that a background source is simultaneously microlensed by 
n-pointlike masses. 

Basically, for an ensemble of compact objects as microlenses which are often assumed to be on a single 
lens plane, the total magnification of a pointlike source at the position /3 is the sum of the magnification /Xj 
of each micro image at position di on the lens plane 



However, this straightforward method cannot be efficiently employed for a computation of the total mag- 
nification. In practice, one is unable to find analytically all the micro images when the number of lenses is 
very large, especially for extended sources. Many numerical techniques have been developed to deal with the 
problem of large number of microlenses including Monte-Carlo simulation (Young, 1981), the ray-shooting 
method (Kayscr, Rcfsdal and StabcU, 1986; Schneider and Weiss, 1987), the Fourier method (Katz, Balbus 
and Paczynski, 1986), the MarkofT method (Deguchi and Watson, 1988), and the parametric representation 
of caustics (Witt, 1990). In particular, the inverse ray-shooting method has been widely used in recent 
years in microlensing simulations: Light rays are traced backwards from the observer to the source plane, 
on which the magnification pattern is represented by the intersection of the rays. So, the number density 
of rays is proportional to the magnification. A typical magnification pattern produced by an ensemble of 
N star-like microlenses is shown in Figure 13. If a background source traverses the magnification regions 
due to either the motion of the source itself or the velocity dispersion of the stars associated with the lens 
galaxies, the apparent luminosity of the source would vary as a function of time. To most observers, this 
kind of variability is something like a "noise" . However, it should be noticed that the microlensing-induccd 
variations could be very dramatic sometimes, which may explain the unusual features associated with some 
special objects. 

It was noticed that the violently variable objects 0846-I-51W1 (quasar), AO 0235-1-164 (BL Lac object) 
and PKS 0537-441 (blazar) might be the results of microlensing (Nottale, 1986; Stickel, Fried and Kiihr, 
1988a, b; 1989). It was even speculated that the variability of apparent magnitude of quasars are partially due 
to microlensing rather than their intrinsic physical processes (Peacock, 1986; Kayser, Refsdal and Stabell, 
1986; Schneider and Weiss, 1987). These arguments have recently been strengthened by Hawkins (1993), 
based on the analysis of a complete sample of ^ 300 quasars selected from their variability over 17 years 
(Hawkins and Veron, 1993; hereafter HV). Some typical light curves from their sample are plotted in Figure 
14 for two high redshift {Zg = 2) and two low redshift {zg = 0.2) quasars, respectively. To investigate whether 
these variabilities are intrinsic to quasars or due to gravitational lensing by compact objects along the lines 




-1 



(53) 
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Figure 14: The typical light curves of four variability-selected quasars by HV. The left panels are high 
redshift quasars at 2s ~ 2 and the right panels are low redshift ones at Zs ~ 0.2. 



Figure 15: The number-magnitude relations for the optically selected quasars (open circles) (BSP) and the 
variability-selected ones (open triangles) (HV) 

of sight, an analysis of the time-varying autocorrelation function was made. It turns out that the timescale 
of quasar luminosity variations decreases with increasing redshift. This is inconsistent with the theory that 
the expansion of the Universe should cause observed timescales to increase linearly with (1 + Zg) due to 
time dilation. Therefore, Hawkins concluded that the quasar variabilities cannot be intrinsic to quasars 
themselves, and gravitational lensing is the most possible cause, indicative of the existence of a large number 
of compact objects up to 10% of the critical mass density of the Universe. Nevertheless, this claim should 
be taken very cautiously in the sense that the observed feature of quasar variability increasing with their 
redshift can also be interpreted as the result of cosmic evolution or various observational limitations, e.g. 
the finite duration of the monitoring campaign, the finite photometric sensitivity (Alexander, 1995) and the 
observing wavelenght dependence arising from the accretion disk model of quasar (Baganoff and Malkan, 
1995). Another interesting issue is the number deficit in HV, as compared with the optically selected quasars 
(Boyle, Shanks and Peterson, 1988; hereafter BSP) (see Figure 15). It remains worth investigating whether 
the quasar number discrepancy in HV is related to microlensing or to the observational methods. Finally, 
the most important aspect of studying cosmological microlensing is to set constraints on the fraction of 
compact objects (f^c) in the matter density of the Universe by analyzing complete samples of variability- 
selected sources, which is quite similar to the purpose of the ongoing MACHO experiments in our Galaxy, 
although one cannot separate the microlensing-induced variability from the variability intrinsic to sources. 
For instances, using the HV sample, Schneider (1993) obtained an upper limits of fic < 0.1 for compact 
objects with masses ranging from 1O~^M0 to 3 x 1O~^M0. Dalcanton et al. (1994) have recently found 
Vtc < 0.1 in the mass range O.OIMq - 2QMq, fl^ < 0.2 for O.OOIM© - 60Mq and < 1 for O.OOIM© - 
3OOM0 by comparing the distributions of the AGN and quasar equivalent widths of emission lines at low and 
high redshifts. It is expected that observations of cosmological microlensing can set more stringent limits on 
Oc in the next few years. 



2.4 Inhomogeneous Universe 

Astrophysical observations indicate that the Universe tends to be locally inhomogeneous on scales less 
than ^ 100 Mpc. Since the early 1960's there have been many studies about the influence of matter 
inhomogeneities on the propagation of light rays from distant sources, especially on the magnitude-redshift 
(m ~ z) relation (Zel'dovich, 1964; Bertotti, 1966; Gunn, 1967; Kantowski, 1969; Dyer and Roeder, 1972; 
1973; Canizares, 1982; Nottale, 1982a,b; 1983; Vietri and Ostriker, 1983; Schneider and Weiss, 1988a,b; 
Isaacson and Canizares; 1989; Wu, 1990b; 1992b; Kantowski, Vaughan and Branch, 1995). Many authors 
addressed the question if the classical m ~ 2; relation (Mattig, 1958) should be modified, i.e., if the m ~ 2; 
relation in a locally inhomogeneous Universe differs from that in a standard Friedmann-Lemaitre Universe. 

In an inhomogeneous Universe the propagation of a bundle of light rays is controlled by two different 
effects: A beam of light traveling outside the mass clump would diverge faster because the matter density in 
such a region is lower than the mean matter density of the Friedmann-Lemaitre Universe, whereas a beam 
of light passing near the clump would be sheared by the gravity of the clump, leading to the convergence of 
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light rays. Flux conservation requires that the divergence due to the absence of matter inside the beam be 
balanced on average by the convergence due to the gravitational effect of the clump, so that the luminosity 
distances in both the inhomogeneous Universe and the Priedmann-Lemaitre Universe remain statistically 
equal if the size of the inhomogeneities is sufficiently small (typically, the size of galaxy) (Weinberg, 1976). 

Suppose that the Universe is uniformly filled by both the intergalactic medium of density of dp and matter 
clumps of density of (1 — a)p so that the mean mass density is the same as that in the Priedmann-Lemaitre 
Universe, a denotes the fraction of the total mass density that is intergalactic. a = 1 corresponds to a 
completely homogeneous Universe, i.e., the Priedmann-Lemaitre model, and a = describes a completely 
inhomogeneous Universe in which all the matter is concentrated into clumps. Por a beam of light rays with 
vertex at the observer (z = 0) propagating far away from any matter clumps, the shearing effect of the beam 
can be neglected and the propagation of light is determined by the optical scalar equation (Sachs, 1961) 

in which we use the angular diameter distance D (in units of c/Hq) as the variable (Dyer and Roeder, 1973), 



is the cosmological mass density parameter and 



O = (55) 



is the critical mass density of the Universe. The initial conditions of eq.(54) can be conveniently chosen to 
be the values of angular diameter distance and the local expansion of the Universe at the observer: 



£'U=o =0; 
4^L=n = l. 



dz 1 2=0 



Under these initial conditions the solution to eq.(54) with a = 1 gives the well-known angular diameter 
distance in the Priedmann-Lemaitre Universe 



- _ 2[zn + {n-2){-i + VT+m)] 



In the case of a = eq.(54) reduces to (Dyer and Roeder, 1972) 



D= / , (59) 



In particular, the solution to eq.(54) for a flat Universe of f2 = 1 is 



2 



D = ^{l + zf^-''^l'' l-(l + z)-'^/" (60) 



in which p = \/25 — 24q. Seitz and Schneider (1994) have recently shown that a general solution to eq.(54) 
can be found for any values of Q and a by transforming eq.(54) into the Legendre differential equation. 

Luminosity distances in the Priedmann-Lemaitre Universe and in the inhomogeneous Universe can be 
denoted by Dlo and Dl, respectively, and relate with angular diameter distances by multiplying a factor 
of (1 + z^. Their difi'erence reflects the divergence of the light propagation in the Universe, whilst the 
convergence can be described by gravitational lensing effect when a beam of light rays passes near the 
clumps. Hence, the parameter 

Am = 5 log ^r-^^ — 2.5 log (61) 



D 



LO 
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indicates a deviation of the actually observed apparent magnitude of a source from its theoretically expected 
value in a completely homogeneous Universe. Note that Am may have relatively large variations, i.e., a 
distant source may dramatically change its apparent luminosity. This arises because the magnification fj, can 
vary in principle from unity to infinity, depending on the distance of the hne of sight to the background source 
from the mass clumps. However, energy conservation requires that luminosity distances in a homogeneous 
Universe and a clumpy Universe should be statistically equal, which reads (Ehlers and Schneider, 1986) 

Thus, the mean magnification correction is 

(Am) =2.5[log(M)-(logM)]. (63) 

The probability that a source at Zs is gravitationally magnified by a factor of /x due to pointlike lenses 
within {zd, Zd + dzd) is 

dp = rid {t^S^DI) {drprop/dzd)dzd (64) 

where is the number density of the lenses, dvp^op, the differential proper distance around the lens at Zd, 
and ■kO'^D'^ is the lensing cross-section which is given by the lensing equation eq.(24): 

Assuming a uniform distribution of pointlike lenses in the Universe, i.e. = (1 + Zd)^ndo, and defining the 
mass density parameter of the lenses as 



UdoM SnGAdndo 



(l_a)0^^= (66) 



we have the total probability for a source at Zg to be magnified by a factor of /x due to foreground compact 
objects 

Pi(M)=30(l-a)^f^^-l) r^4^dzd. (67) 




Fortunately, this expression can be separated into two parts 

Pi(m) = hin) T (68) 

where 

/i(m) = 2(^^-1); (69) 



^0/1 ~\^o f^" DdDds l + Zd , 



T is the optical depth to gravitational lensing in terms of the definition of eq.(45), which describes the total 
number of lenses enclosed within the Einstein ring along the line of sight to the source. Figure 16 shows the 
"maximum" optical depth contributed by pointlike lenses that compose all the matter of the Universe, i.e., 
a = 0. Note that eq.(70) is valid for computation of the optical depth to gravitational lensing by various 
compact objects with different masses. Meanwhile, r is an indicator of the significance of multiple lenses. 
In the case of a = and = 1 (Figure 16), a source with Zg > 2> may be lensed by more than one lensing 
object (r ^ 1). Therefore, multiple lenses need to be taken into account for the high-redshift sources if a 
relatively large amount of matter of the Universe is concentrated into compact objects. 

Multiple lenses fall into two classes: geometrically-thin multiple lenses located on a single lens plane, 
as was discussed in the above subsection, and spatially discrete multiple lenses which are distributed on 
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Figure 16: The optical depth to gravitational lensing by the pointlike lenses in a completely inhomogeneous 
Universe ((5 = and f2 = 1). 

different lens planes along the line of sight. For the multiple lens plane deflections, the gravitational lensing 
equation can be formally written as (Blandford and Narayan, 1986; Schneider, Ehlers and Falco, 1992) 

D ^ 



i=l 

i-1 



where ij is the position of the source in source plane, Q;j(^j) is the deflection angle of the light ray with 
the impact distance by the deflectors in the i-th lens plane, £>y denotes the angular diameter distance 
from the i-th lens plane to the j-th lens plane and Di, from the observer to the i-th lens plane (Note that 
j = s refers to the source). The total magnific;atioii of the source luminosity is finally given by the inverse 
of the determinant of the magnification matrix which relates to the Jacobian matrices of the mapping from 
the {i — l)-th lens plane to the i-th lens plane (1 < i < N). Both numerical techniques (e.g. Schneider 
and Weiss, 1988a, b) and analytical methods (Seitz and Schneider, 1992, 1994) have been employed in the 
determination of the magnification factor as well as the magnification probability of a distant source by the 
deflectors in multiple lens planes, although the procedures turn to be relatively complex. 

In the practical computation of magnification probability by miiltiple lenses, some approximations are 
often employed in order to avoid the above complexity. The fact that the probability Pi {/i) by a single lens 
factorizes into two independent parts in eq.(68) leads to the speculation that the probability P„(/u) for a 
distant source by n lenses be written as a product (Wu, 1990a) 

Pnifi) = /n(M) 9n{r), (73) 

i.e.. the /x variable separates from the r variable. gn{T) is a function that describes the probability of a 
background source being lensed by n foreground objects. Obviously, gn{T) follows a Poisson distribution 

9n{r) = J e-\ (74) 
Therefore, the magnification distribution function can be obtained by summing up Pnilj) over all n 

^(z^) = E ^^^^ 

Unfortunately, no exact expressions have been established for P>2(^). An oversimple assumption is that the 
total magnification /i is the product of the individual magnification /i^ so that P(/i) is the convolution of 
Pi{fi) of each lens (Canizares, 1982; Vietri and Ostriker, 1983; Peacock, 1986; Isaacson and Canizares, 1989; 
Schneider, 1993; Pei, 1993). Nevertheless, a justification for this assumption is rather hard, and it seems 
that this multiplication method may be questionable (Wu, 1990a). 

Finally, in a single lensing approximation the mean correction of the apparent magnitude [eq.(63)] to the 
Mattig's relation is (Wu, 1990b,1992b) 

(Am) = 1.16t -t- 0{t^) « 0.58(1 - a)in/2)z'^ + • • •, (76) 

in comparison with the correction found by Kantowski (1969) and Dyer and Roeder (1974) based on the 
Swiss-cheese model for the inhomogeneities of the Universe: 

(Am) w 1.086(1 - a){n/2)z'^ + ■■■. (77) 
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Figure 17: Gravitational lens system QSO 0957+561A,B. [from E. Schild (1991)] 



Figure 18: VLBI observation of the gravitational lens system QSO 0957+561A,B. [from Garrett, et al. 
(1994)] 

These formulae can be used to statisticahy estimate the modification to the classical magnitude-redshift 
relation due to the matter inhomogeneities, and the quantity depends sharply on the matter content of 
compact objects (1 — a)0 in the Universe. 



3 GALAXIES AND MULTIPLE IMAGES 
3.1 Multiply- imaged quasars 

The most significant feature of gravitational lensing is the multiple imaging of the lensed source. The 
first gravitational lens system was discovered in 1979 (Walsh, Carswell and Weymann, 1979) during the 
identification of the optical counterpart of a radio source. This famous system 0957+561 consists of two 
quasar images (A and B) at the same redshift of Zg = 1.41 and with separation of 6".l. A galaxy at Zd = 0.36 
was soon detected at the position near the image B (Adams and Boroson, 1979; Yoimg et al., 1980; Stockton, 
1980), which is believed to be the main defiector with mass of ~ IO^^Mq (see Figure 17). VLBI observation 
of 0957+561A,B (e.g. Gorenstein et al., 1988; Garret et al., 1994) strongly confirmed the lensing origin of 
QSO 0957+561A,B by revealing the same radio morphology (one compact core and three jets) of the two 
images (Figure 18). 

The criteria for the determination of a gravitationally lensed quasar system are (l)multiple images, 
(2)similar spectra and the same redshifts of the images and (3)the detection of intervening galaxies as 
lenses. A list of the accepted and the proposed lensed quasar systems compiled by Surdej and Soucail in 
1993 is updated in Table 3. The proposed rather than confirmed cases arise mainly from the fact that 
no corresponding deflectors have been found to be responsible for the multiple images. The absence of the 
lensing galaxies in the direct imaging centered on some of the multiple quasars even with HST is a well-known 
puzzle in gravitational lensing. It remains unclear today if one should really reject these lensing candidates 
in which the "dark" lensing galaxies are apparently missing. The recent detections of a very faint galaxy 
{B = 25.0) and the faint cluster of galaxies at z > 1 (Mellier et al., 1994) associated with the double quasar 
2345+007 (Fischer et al., 1994) may be promising for the future searches of the deflectors in other lensing 
candidate systems. 



Table 3: Gravitationally-Lensed Multiple Quasars 
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Figure 19: 15 GHz image of the radio Einstein ring MG1131+0456 [from J. Hewitt et al. (1988)] 



3.2 Multiply- imaged radio sources and radio rings 

Radio observation turns out to be an efficient way of finding multiple images due to its consistently high 
dynamical range and resolution of the maps. Indeed, the VLBI observation of double quasar 0957+561A,B 
(Figure 18) played an important rok; in the confirmation of their Icnsing origin (Gorenstcin ct al., 1988). 
Today, with the completeness of a large lens survey at the Very Large Array, the Cosmic Lens All-Sky Survey 
would provide a large number of lensing candidates for both modeling of the lensing systems and statistical 
study. The recent discoveries of a double lens images 1600+434 (Jackson c;t al., 1995) and a quadruple lens 
system 1608+656 (Myers et al., 1995; Snellen et al., 1995) have marked the success of this survey. 

Another success of radio observation in gravitational lensing is the detection of the Einstein ring. As was 
illustrated in Figure 2, when the background source lies in a position behind the foreground lensing object, 
the multiple images merge into a ring-like image, i.e. the Einstein ring. The first radio ring (Figure 19) was 
discovered by Hewitt et al. in 1988, which is the image of a radio lobe at Zs = 1.13 lying perfectly behind a 
foreground galaxy at Zd — 0.85. Six more radio rings have been so far observed: MG 1654+ 1346 (Langston 
et al, 1989), PKS 1830 - 211 (Rao and Subrahmanyan, 1988; Jauncey et al., 1991), MG 1549 + 3047, MG 
0751 + 2716, B 1938 + 666 (Lehar et al., 1993), and B 0218 + 357 A-B (Patnaik et al., 1993). Besides 
constraining the mass profile in the lensing galaxies from modeling of the radio rings, the measurement of 
time delay in the radio ring would be of great interest for the determination of the Hubble constant (see 
section 3.4). 



3.3 Galetxies as Lenses 

Galaxies are found to be the main deflectors for the ^ 10 confirmed gravitationally-lensed multiple quasars in 
Table 3. Modeling each lens system based on the observational data has reproduced quite well the observed 
multiple images of quasars. Actually, this procedure is no more than to solve the lensing equation for different 
gravitational potentials, which has been extensively discussed in section 1.2. The most important issue in the 
study of multiply-imaged quasars, however, is the statistical properties, which raises the question if there 
are enough massive galaxies in the Universe to be responsible for the observed events of gravitationally- 
lensed quasars. Recall that the lensing galaxies have not been detected today in some of the proposed lens 
candidates. Moreover, statistical lensing of galaxies provides the information on how many multiply-imaged 
quasars would be expected to observe over the sky. By comparing the theoretical statistical predictions 
with the observations of multiple quasars, one can also set useful constraints on the mass density of the 
lensing objects in the Universe (e.g., Hewitt et al., 1986; Claeskens et al., 1993; Surdej et al., 1995) and the 
cosmological constant A (e.g.. Turner, 1990; Fukugita et al., 1992; Sasaki and Takahara, 1993; Rix et al., 
1994; Kochanek, 1992;1993a,c;1995). 

Adopting the simplest matter distribution, SIS, for the lensing galaxy, one can write the lensing cross- 
section from eq.(27) to be 



2 



Hence, the probability that a source at Zg is magnified by a factor of greater than /x due to an ensemble of 
galaxies within redshift dzd of Zd is 

dp = F \ — = ^ dzd 7 (79) 

\ D,, J VI + i^zd (M-1)' ^ ^ 

where F = 16TT^no{c/ Hq)^ (ay / c)^ , no is the present comoving number density of the galaxies assumed 
to develop as n = (1 + Zd)^nQ, Dd, Dg and Dds are the angular diameter distances in units of {c/Hq), 
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Figure 20: The differential optical depth to gravitational lensing by an ensemble of SIS galaxies. The 
cosmological model is chosen to be f2 = 1 and a = 1. F = 167r^no(c/iJo)^(ft,/c)'^. (cf. Turner, Ostriker and 
Gott, 1984) 



corresponding to Dj^, Ds and Dj^g, respectively. In a way similar to the lensing probability by pointlike 
masses, eq.(79) can be separated into two parts by utilizing the optical depth r 



1 



* = 7—^- (80) 



Figure 20 shows the differential optical depth dr to gravitational lensing for a source at different redshifts 
Zg = 1,2,3. The significance of the differential optical depth is that it provides a clear view of the most 
probable lens position for various sources. It then turns out that for quasars at their typical position of w 2 
the lensing galaxies locate most likely at Zd ~ 0.5. The total optical depth r can be obtained analytically 
for some special cosmological models (Turner, Ostriker and Gott, 1984). For example, in the case of O = 1 
and a = 1 the total optical depth is simply 

4F [{l + Zg)'/^-ir 
" = l5 (l + ..)3/2 • (''^ 

For a quasar a,t Zg = 2 this reads 

r = 2 X 10"^(F/0.1) (82) 
An extensive analysis of various types of galaxies by Fukugita and Turner (1991) gives 

F = 0.019 ± 0.008 for E galaxies 

F = 0.021 ± 0.009 for SO galaxies (83) 
F = 0.007 ±0.003 for S galaxies 

based on a morphological composition E:S0:S=12: 19:69. Thus, for all the galaxies as lenses, F = 0.047. This 

results in a total optical depth of 0.001 for a quasar at Zg = 2, i.e., about 1/1000 quasars at Zs ^ 2 would 
be found to be significantly lensed by foreground galaxies. Note that the definition of the optical depth 
utilizes a cross-section of n6^, so that the total optical depth in SIS corresponds to the total probability of 
a background source being magnified by a factor of /i > 2. Using other models (ISC, KING, r^^^, etc.) for 
the matter distribution of galaxies has yielded statistical properties which explain very well the observed 
frequency of lensed quasars, the observed distributions of image separations and of apparent magnitudes 
(Dyer, 1984; Hinshaw and Krauss, 1987; Kochanek and Blandford, 1987; Wu, 1989b; Mao, 1991; Fukugita 
and Turner, 1991; Kochanek, 1993a,b,c;1995; etc.). 

It is worth noticing that the Hubble Space Telescope (HST) Snapshot Survey (Bahcall et al., 1992; 
Maoz et al., 1992; Maoz et al., 1993a, b; Falco, 1993) provides a sample of 502 luminous and high redshift 
quasars, among which a search for gravitationally lensed events has been made using the HST Planetary 
Camera. This sample is of great significance for testing the theoretical model of gravitational lensing which 
predicts that there should be many gravitationally lensed quasars having ^ 1 arcsecond (Turner, Ostriker 
and Gott, 1984) and even sub-arcsccond image separations (Fukugita and Turner, 1991), whilst the HST 
Snapshot Survey is capable of observing these small separated cases. One new candidate Q1208+1011 was 
found in the HST Snapshot Survey. Together with the previously known cases in the sample, the observed 
frequency of lensing is estimated to be between 3 and 6 out of 502 quasars. The theoretically expected 
frequency of lensing can be obtained by T{Zg)B{m, z), where B{m, z) is the factor by which lensed quasars 
are over-represented among quasars of a given magnitude m because fainter quasars have been magnified 
to that magnitude, and its value is determined by the quasar luminosity function (Bahcall et al. 1992). As 
a consequence, this indeed results in a frequency compatible with the result of the HST Snapshot Survey 
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but does not meet the prediction by a cosmological constant dominated Universe (Maoz et al., 1993b). The 
newly completed lens quasar surveys have even limited Aq = A/3Hq to Aq < 0.66 (Kochanek, 1995). 



3.4 Determination of Hq 

The present status of the uncertainty of determination of the Hubble constant Hq by a factor of about two 
between 40 km/s/Mpc and 80 km/s/Mpc is unfortunate, which has caused the cosmic distance dispute for 
decades. The main problem arises from the disagreement of the "standard candles" used as the indicator of 
the absolute distance (see Fukugita, Hogan and Peebles, 1993 for a recent review). Indeed, it is very unlikely 
that the debate on Hq would be settled in the next few years if the measurements are still based on the 
conventional "standard candle" methods. However, the recent progress of determination of Hq using other 
techniques that arc independent of the usual distance-ladder arguments may hopefully help to settle down 
the debate. The time delay of the multiply-imaged quasars is one of these methods which are promising for 
the measurement of Hq. 

The time that the photon traverses a proper distance of Arprop in the Universe is simply 

^ Arp^p ^ 1 Az ^ 1 

Therefore, wo can obtain the Hubble constant by measuring the time difference At for a given Q. The 
different optical paths to earth between the multiple images of a lensed quasar offer then a possibility of 
observing At (the time delay) if the quasar has intrinsic luminosity variability. 

Actually, even before the discovery of the gravitationally-doubled quasars it was realized that the Hubble 
constant Hq could be measured from the time delay of two images of a single background source (Refsdal, 
1964a,b, 1966). The expression eq.(84) can now be generally written as 

HoAt = T{n,zd,zs) h.ns{eA,B,a) (85) 

where T is called the cosmological correction function which is only dependent on the cosmological parame- 
ters, and /lens is the lens model function which is given by the matter distribution of the lens. 9 and a are 
the positions of the images and the deflection of light, respectively. It has been shown (Kayser, 1986) that 
such a separation of the cosmological function T from the lens model function / is indeed possible. 

There are essentially three approaches developed to compute the time delay, namely, the wavefront 
method (Kayser and Refsdal, 1983); (2)the integration method (Cooke and Kantowski, 1975) and (3)the 
scalar formulation (Schneider, 1985). In the approximation of weak gravitational field (^), the traveling 
time f of a photon is (Cooke and Kantowski, 1975; Borgeest, 1983) 



ct 



= j ds - ^ j (t)ds. (86) 



The first term gives the length of the light path, and the second is the relativistic time dilatation due 
to the gravitational potential 4> of the dc;flcctor. The integrals are performed along the photon orbit ds. 
Correspondingly, the propagation delay for light-rays from the double images of a lensed source can be split 
into two components: (l)the geometrical delay Atg that is caused by the difference of light paths between 
images and (2)the potential delay Atp that is induced by the difference of the gravitational field of the 
intervening lensing object at the image positions. A straightforward computation gives the two components 
(Cooke and Kantowski, 1975; Borgeest, 1983; Borgeest and Refsdal, 1984): 

cAt, = (1 + ..)iP.^M±^ • {6, - e.) (87) 



cAtp = (1 + Zd)Dd / a(6/) • dd. 
J Ox 
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Figure 21: Optical light curves of the double images of QSO0957+561. (From Vanderriest et al., 1989) 



Figure 22: The updated optical light curves of the double images of QSO0957+561 A (upper) and B (lower) 
with the B data retarded by 1.1 years and offset by 0.2 mag. The 832 actual observations of the image A 
are shown at the bottom. (From Schild and Thomson, 1995) 

Finally, extracting Hq from Dgi and combining it with At — Atg + Atp yields the form of eq.(85). 

A detailed lens model has been developed for the well-known lens system QSO0957+561 in order to 
determine the Hubble constant Hq from the gradually accumulated data of the light curves of the double 
quasar images A,B (Young et al., 1980; Dyer and Roeder, 1980; Greenfield, Roberts and Burke, 1985; Falco, 
Gorenstein and Shapiro, 1985;1991; Gorenstcin, Falco and Shapiro, 1988; etc.). A lensing model that is com- 
posed of three matter components has been found to reproduce the known properties of QSO0957+561A,B 
quite well: (l)thc bright galaxy (Gl) described by a King profile ((Jt,, Tc); (2)a compact nucleus with mass 
of Mc ^ lO^^h^^ Mq and (3)the cluster characterized by the surface mass density of a smoothly distributed 
mass screen (Gorenstein, Falco and Shapiro, 1988; Falco, Gorenstein and Shapiro, 1991). This lens model 
leads to a relatively simple expression of the Hubble constant Hq from eqs.(85), (87) and (88) (see also 
Roberts et al., 1991) 

The error estimate includes measurement error between the VLBI, VLA and optical observations, the un- 
known values of Q and of the dumpiness of the Universe, the non-uniqueness of the cluster model, and errors 
in the detailed model of Gl. The result depends weakly on the cosmological density parameter fl. It appears 
that two free parameters, and At, control the actual evaluation of Hq. Motivated by the significance 
of determination of Hq from the time delay, Rhee (1991) obtained the line-of-sight velocity dispersion of 
303 ± 50 km/s for the bright galaxy Gl in the QSO0957-I-561 lens system, leaving the final work of finding 
Hq to the measurement of the time delay between the double images. However, it should be mentioned that 
the simple lens model from which eq.(89) was derived is by no means unique, and more complicated mass 
distributions are not only possible but also actually well motivated (Bernstein, Tyson and Kochanek, 1993). 

Vanderriest et al. (1989) undertook a 8-years optical photometric monitoring of the double quasar 
0957-|-561A,B from 1980 to 1987, which contains totally 131 observations. The light curves of the image A 
and B are shown in Figure 21. A significant decrease in brightness around Julian Days 2445700 in A and 
Julian Days 2446100 in B is clearly seen, indicative of a time delay ~ 400 days. The original analysis of 
Vanderriest et al. (1989) from the cross-correlation function for the two light curves gives At = 415 days, 
while Press, Rybicki and Hewitt (1992) reached a value of At = (537 ± 11) days based on a newly developed 
mathematical methodology for the same data. Using additional optical data of 3.5 years coverage, Schild 
(1990) obtained a value of At = 404 days, consistent with the result of Vanderriest et al. (1990). The 
updated optical light curves to 1994 July (Figure 22) seems also to support At « 1.1 years (Schild and 
Thomson, 1995). Radio monitoring of QSO0957-|-561A,B with the VLA was reported later for a 10 - 11 
years coverage (Roberts et al., 1991; Lehar et al., 1992). The fiux curves of the two images are shown in 
Figure 23. In the absence of specific features in the two curves, an analysis of the cross-correlation function 
of two signals suggests a time delay of 513 ± 40 days. The very recent measurement of the time delay using 
the hybrid maps of QSO0957+561A,B with VLBI spanning the 6-year interval (1987-1993) yields At - 1 
year (Campbell et al., 1995), while a reanalysis of the Lehar et al. (1992) observation with a refined method 
has found that their radio data are compatible with the result of 1 year obtained from the optical data 
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Figure 23: Radio flux curves of the A (uppercase) and B (lowercase) images of the double QSO0957+561. 
The letters denote the VLA configurations (P for partial configurations). (From Lehar et al., 1992) 

(Pelt et al. 1995). Moreover, the updated VLA light curves of QSO0957+561A,B for 16 years coverage show 

a time delay of At = 455 ± 40 days (Haarsma ct al., 1995). Apparently, the present results of measurement 
of the time delay in the double images QSO0957+561A,B are controversial, and the acceptable value ranges 
from 404 days to 537 days. 

Adopting the cosmological density parameter of 17 = 1 and the velocity dispersion of cr„ = 303 ±50 km/s 
for the Gl in QSO0957+561A,B system, one finds from cq.(89) 

u f 48+7*^' 1 , / /TVT f At = 415 ± 20 days (optical) ) 

= I 39i I ^"^/^/^P^ { At = 513 ± 40 days (radio) | " ^''^ 

Though these values still contain large uncertainties, they seem to support a low value of the Hubble con- 
stant. Further observations will be needed to find a reliable value of time delay in QSO0957+561A,B as well 
as to establish a reliable lensing model in order to precisely determine the value of Hq. The new observation 
of luminosity variations in the quadruple-lens system B1422-I-231 (Hjorth et al, 1995) and the recent detec- 
tion of time delay in the Einstein ring B 0218-1-367 (Corbett et al., 1995) and PKS 1830-211 (van Ommen 
and Preston, 1995) would be also promising for the measurement of Hq. 



3.5 Quasar-galaxy associations 

One of the important consequences of gravitational lensing, as first realized by Gott and Gimn (1974) even 
before the discovery of the first gravitationally-imaged quasar 0957+56lA,B, is that the surface number 
density of quasars near foreground galaxies would be enhanced (denoted by the quasar enhancement factor 
qq) becaiisc the distant quasars lying behind galaxies would be magnified by the lensing effect of the galaxies 
and then enter into the detection limit (see also Canizares, 1981; Vietri and Ostriker, 1983; Schneider, 
1986;1987a,b; Kovner, 1989; etc.). Equivalently, an overdensity of foreground galaxies around high-redshift 
quasars would also exist (described by the galaxy enhancement factor qo) (Schneider, 1989). The statistical 
evidence on such quasar-galaxy associations was firstly found by Tyson (1986) and later reported by Webster 
et al. (1988). They all claimed a significant enhancement of galaxy surface density in the vicinity of 
distant quasars. Since then, the observational evidences for quasar-galaxy associations have been cumulated 
(Narayan, 1992). Table 4 summarizes the present status on the optically-selected quasar-galaxy associations, 
including two negative results. Some suggestions have been made to improve the confidence of the different 
results such as choosing the same objects, cross-calibrating the different observing techniques, using the same 
criteria, etc. Yet, large samples will be needed to further confirm the existence of quasar-galaxy associations. 

The first effect of gravitational lensing is its magnification (/i), which enhances the apparent brightness 
of background sources by an amount of 2.5 log /x in magnitude, leading to an increase of the surface number 
density (a) of background sources by picking up the faint sources: a{6) N{< m + 2.5logfi)/So{6), where 
50(6*) is the observed area at a distance 6 from the deflector. The second effect is the area distortion 
So {0) + ASq (0) , which arises from the light bending around the deflector (see Figure 24) . This reduces the 
number counts by losing the sources within the dashed-line regions: a{9) ~ N{< m)/{So{6) + ASo{6)). As 
a whole, the surface number density can be written as 

N{<m + 2.5\og^i) 
^^'^ = S,[0) + ASo[0) ■ (^1) 

Defining the enhancement factor qQ{0) as the ratio of the disturbed surface number density a{0) to the 
undisturbed one iJo(^) = A''(< m)/So{0) and noticing that /x(^) = [5o(^)-|- AS'o(^)]/6'o(^), one has (Narayan, 
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Table 4: Foreground galaxy enhancement qc 



authors 


QSO 


selections 


9 range(") 


galaxy(i?) 


QG 


Crampton 


101 


V < 18.5 

z > 1.5 


0-6 


~ 23 


1.4 ±0.5 


Kedziora— 


181 


V < 18.5 


6-90 


~ 21.5 


~ 1 


Chudczer 




z > 0.65 








Magain 


153 


V = 17.4 
(z) = 2.3 


0-3 


~ 21 


~ 2.8 


Thomas 


64 


V < 18.5 

1 < z < 2.5 


0-10 


~ 22 


~ 1.7 


Van Drom 


136 


V = 17 A 
(z) = 2.3 


3-13.7 


~ 23 


~ 1.46 


Webster 


68 


V <18 
0.7<z<2.3 


3-10 


~ 22 


~ 2 


Yee 


94 


V <19 
z > 1.5 


2-6 
2-10 
2-15 


~ 22.5 


1.0 ±0.3 
1.0 ±0.2 
0.9 ±0.1 



Figure 24: Scheme of two effects of gravitational lensing on the number counts. Magnification effect enhances 
the apparent magnitude of the background sources, which helps to pick up the fainter sources and then leads 
to an increase of the total number of the sources in a flux-limited sample. The area distortion effect due to 
the light deflection reduces the source volume (the dashed-line region), resulting in a decrease of the total 
number counts. 
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Figure 25: The local quasar enhanccnicnt qq against quasar limiting magnitude B and local lensing mag- 
nification fj, for the BSP and the HV counts. The solid lines correspond to the results within the BSP and 
HV survey limit of -B = 21, i.e., both B and B + 2.5 log /i in the computation of qq are confined to B < 21, 
and the dotted lines are the extrapolated results by employing A^(< B) beyond B = 21. Note that in the 
HV counts, qg > 1, providing always the positive associations, while in the BSP counts qg may be smaller 
than unity, leading to the "negative" associations. 



1989) 

^^^^^ a^Rh WY ^ ^ 

It appears that the "local" enhancement parameter qg at 9 depends on two factors: the intrinsic number- 
magnitude relation of the background sources and the local lensing magnification around 9. The average 
quasar enhancement qg(9\,62) over an angular distance of (^1,^2) from the foreground deflector is simply 

qg{9i,e2)= (g2_g2) • (93) 

Two kinds of quasar number-magnitude relations have been thus far suggested from, respectively, the 
BSP survey and the HV survey. The significant difference in these two relations is that the slope of the BSP 
counts changes at B w 19.15 from 0.86 to 0.28 while there is no such a turnover in the range of B < 21 in 
the HV counts (see Figure 15). This probably arises from the different selection methods used in the two 
surveys. The BSP cumulative counts can be fitted by (Narayan, 1989) 

iV(< B) = 4.66 X lOO-s^^^-i^ i'"^), B< 19.15; , . 

iV(< B) = -10.95 + 15.61 X ioO-28(B-i9.i5). q > 1915. 

This relation is valid for z < 2.2 and i? < 21. Nevertheless, the subsequent observation (Boyle, Jones and 
Shanks, 1991) indicates that the above relation holds true also to B < 22. The HV cumulative counts can 
be fitted by (Wu, 1994a) 

7V(<B)=6.25xl0°-5i(^-i9i5). (95) 

The local enhancements qg are computed for these two kinds of quasar number counts and plotted 
against the local magnification /n in Figure 25. One should pay a special attention to the case where 
B + 2.5 log /i is larger than the limit of validity of the quasar number count relation. BSP and HV surveys 
were both restricted within B < 21. Therefore, the number-magnitude relation A^(< B + 2.5 log /i) fails 
when B + 2.5 log > 21, which occurs for a sufficiently large fi. Strictly speaking, one cannot calculate the 
enhancement qg beyond the survey limit, and the extrapolation of the solid lines in Figure 25 requires the 
knowledge of fainter quasar counts. Furthermore, the application of eqs.(94) and (95) for the evaluation of 
qg in eq.(92) has presumed that the observed iV(< B) remains the same as the intrinsic counts, i.e., quasar 
counts have not been contaminated significantly by lensing. 

A power-law number-magnitude relation with index of a, N{< m) ^ 10"™, would lead to qg = /i^ "''"^^, 
independent of the limiting magnitude m. The HV counts have log qg/ log fi — 0.3, which then cannot 
provide a large enhancement for a moderate magnification. Conversely, the two power-laws of BSP counts 
give rise to a relatively wide range of qg, depending on both the limiting magnitude and the magnification. 
Bright quasars {B < 18) appear to be relatively strongly associated with foreground deflectors, with a 
maximum enhancement at dqg/dix = 0. On the other hand, faint quasars {B > 19.15) exhibit a "negative" 
association with the foreground deflectors, i.e., fewer quasars would be found near the foreground galaxies 
than in the rest of the sky. 
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Figure 26: The average quasar enhancement qq over different search areas around a foreground galaxy with 
(Jv = 215 km/s. 

Adopting a SIS model for the matter distribution of a lensing galaxy, we can express the local magnifi- 
cation at an angular distance 6 from the center of the galaxy as [eq.(27)] 



Furthermore, wc assume the distance parameter D^s/ to be very close to unity, which approximately holds 
true when the foreground galaxies are at relatively low redshift while quasars are at high redshift in order 
to guarantee that they are not physically associated systems in the searches for quasar-galaxy associations. 
In this case, the Einstein radius reads 



Figure 26 shows the average enhancement qq versus the limiting magnitudes and the search ranges around 
a typical galaxy of ct = 215 km/s, the average of the E/SO galaxy velocity dispersions (Kochanek, 1993c), 
provided that the extrapolation of both BSP and HV counts to the faint magnitude {B > 21) is possible. In 
HV counts qq depends only on the search areas and the resulting amplitude turns to be too small to explain 
the reported enhancements of as large as 2 listed in Table 4. Both positive and "negative" associations 
are provided by BSP counts, separated in the range of 19 < B < 20. Other important conclusions are: 
(l)Positive associations between foreground galaxies and background quasars would be found when one 
chooses the limiting magnitude of the quasar sample to be brighter than B « 19. (2)When the faint 
quasars {B > 19.5) are involved, one would expect to detect null or "negative" associations. This scenario 
of existence of positive/negative quasar-galaxy associations can explain the observed results (Table 4) quite 
well (Wu, 1994a). In fact, the search for "negative" associations between quasars and galaxies at faint 
magnitude can be used as the ultimate test for whether or not the quasar-galaxy associations stem from the 
effect of gravitational lensing. There are no other mechanism known thus far that can result in the negative 
associations, i.e., the surface number density of distant quasars around foreground galaxies is smaller than 
the mean quasar density in the rest of the sky. 

For foreground galaxies with different luminosities which follow the Schechter function (j){L)dL = (j>*{L/L^) 
exp(— L/L«)(i(L/L*), we can find the average quasar enhancement over a search range of (0i,02) from 



where i and ji represent, respectively, the i-th type and composition of galaxies, and the minimum galaxy 

luminosity Lmin is related to the galaxy limiting magnitude in the observation of quasar-galaxy associations. 
The galaxy luminosity L can be converted into the velocity dispersion ay through the empirical formula such 
as the Faber-Jackson relation for early-type galaxies (E/SO) L/L^ = (ay/a^)^ or the Tully-Fisher relation 
for spiral galaxies (S) L/L^, = ((Tt,/(T*)^'^. The parameters L* (or cr*). (f>* and z/ have been observationally 
determined for different galaxies (see Fukugita and Turner, 1991). Finally, the spatial distribution of galaxies 
needs to be taken into account. Assuming a constant comoving number density of galaxies, we have 




(97) 




(98) 



{qQ{0i,02)) = 



/;° AnPlil + z^f jqN) dr, 



prop,Zd 



prop,Zd 



(99) 
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Figure 27: Giant luminous arcs in Abell 370 (a) (Figure courtesy of G. Soucail) and CI 2244-02 (b) (Figure 
courtesy of F. Hammer) 

Table 5: Arclike Images 



in which drprop,za = {c/ HQ)dzd/[{l + \/T+Tlzd\. This theoretical expectation can be used straightfor- 
wardly to compare with the observations (Wu, Zhu & Fang, 1995). 

It should be realized that the quasar-galaxy associations are the statistical results and thereby, statistical 
lensing should be involved. An exact treatment of this question is to convolve the magnification probability 
P(/z) by foreground galaxies with the intrinsic quasar number counts iV(< m) (Schneider, 1989). The 
difficulty is that P(/i) should be artificially truncated at the faint end of magnification in order to perform 
the integration / A'^(< m)dP{n). 



4 CLUSTERS OF GALAXIES AND ARCLIKE IMAGES 
4.1 Giant arcs and arclets 

The arclike image associated with Abell cluster 370 was first detected by Hoag in 1981. However, this blue 
arc had not been recognized to be the imago of a distant galaxy gravitationally lensed by the cluster until a 
few years later when two groups of astronomers independently announced their convincing evidences of the 
existence of this peculiar feature in the Universe (Soucail et al., 1987a; Lynds and Petrosian, 1986), followed 
by Paczynski's (1987) lensing interpretation. Only two giant arcs (A370 and C12244-G2, see Figure 27) were 
known at that time and other explanations remained also possible. The crucial point of interpretation of the 
arcs as gravitationally imaged background sources rather than some peculiar features physically associated 
with the clusters is the measurements of redshift of the two giant arcs. They do show higher rcdshifts than 
those of their associated clusters: The giant arc in Abell 370 has a redshift of 0.725 (Soucail et al., 1987b; 
Miller and Goodrich, 1988), in comparison with the redshift of 0.374 for Abell 370 itself , and the arc in CI 
2244-02 has probably an even higher redshift of 2.237 (McUier et al.. 1991) while its associated cluster is 
only at Zd = 0.336. These measurements have strongly confirmed the lensing origin of arclike images seen in 
the cores of rich galaxy clusters. Up to now giant arcs and arclets have been detected in about 30 clusters of 
galaxies (Table 5) and this number is still increasing dramatically. In particular, the high X-ray luminosity 
clusters in the EMSS sample turn to be the very efficient deflectors of producing arcs and about 14 arcs have 
thus far been seen in a subsample of 41 EMSS clusters with >'2 x 10*^erg/s and za > 0.15 (Hammer et 
al.. 1993; Lc Fcvre et al., 1994; Gioia and Luppino, 1994; Hammer, 1995; Luppino et al., 1995). Note that 
there are several multiple arc systems. 

Wu and Hammer (1993) classified the elongated images using two parameters: axial ratio {L/W, i.e., 
length/ width) and apparent magnitude (B). The "giant" arcs refer to those images whose axial ratios are 
greater than 10, i.e., L/W > 10, and mini-arcs or arclets have L/W < 3. The rest arclike images in between 
(3 < L/W < 10) are called medium arcs. Arc brightness is represented by its B magnitude so that the 
"luminous" arcs have B < 22.5. In the updated list of arclike images of Table 5 there are totally ~ 10 giant 
luminous arcs, while arclets are numerous but usually very faint {B ~ 26). 
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Figure 28: (a)The projected galaxy distribution on the sky without Icnsing and (b) the distorted and 
magnified images of background galaxies by a galaxy cluster (at center) with = 1000 km/s and at 
Zd = 0.25. 



Figure 29: The relative distribution (normalized at LjW = 10) of arcs predicted in Figure 28(b) against 
their axial ratio. 

4.2 Clusters as lenses 

Although most of the arcliko images remain to be spatially unresolved in width today, the colours and spc;c;tra 
of the arcs are compatible with those of local sub-L* spiral galaxies, indicating that they are probable spirals 
at relatively high redshift ~ 1. The foreground clusters often have high X-ray luminosity {L^ > 10*^ 
erg/s) and/or large velocity dispersion (a^ > 1000 km/s), which are then massive enough to act as strong 
lenses for the background galaxies. 

To demonstrate how a rich galaxy cluster at intermediate redshift gravitationally distorts background 
galaxies, a simulation is made based on the current knowledge of dynamical and spatial properties of galaxies 
and rich clusters of galaxies. SIS is adopted for the matter distribution of a galaxy cluster with = 
1000 km/s at Zd = 0.25. The luminosity function established by Broadhurst, Ellis and Shanks (1989) in 
deep redshift survey is employed for the distribution of background spiral galaxies. The luminous area 
of each galaxy is taken to be a circular disk of radius of R which is assumed to follow the relation R = 
i?*(Z//i*)^/^, where i?* is the characteristic radius corresponding to a Z/*-galaxy. In the actual simulation, 
R = io(~i'^ i6~^)/5 kpc (Freeman, 1970). Figure 28(a) shows the unperturbed background galaxies projected 
on the sky in a field of 1' x 1', in which the orientations of the disk galaxies are randomly placed in space 
and the population of the galaxies is presumed not to evolve cosmologically. Moreover, we truncate the 
redshift of the galaxies at = 1.25 duo to the failure of the _ft'-corrcction in the adopted luminosity fimction 
(Broadhurst, Ellis and Shanks, 1989). Figure 28(b) illustrates the same field with a galaxy cluster at the 
center. All the images brighter than B = 23 are selected by taking the magnification effect into account. It 
appears that galaxies in the field have been strongly elongated around the cluster, indicating that distant 
rich clusters of galaxies can indeed act as strong lenses and produce the images of giant arcs and arclets. 

An interesting issue is the total number of mini-arcs and medium arcs with respect to the number of 
giant arcs seen up to the same flux threshold in a galaxy cluster. A statistical distribution of arc number 
appearing in Figure 28(b) against axial ratio L/W is plotted in Figure 29, in which the normalization is 
made at L/W = 10. It turns out that for one observed giant arc, there should be another one medium arc of 
L/W « 6 and two arclets of L/W w 3. The more careful statistical investigations have reached essentially 
similar conclusion (Grossman and Narayan, 1988; Wu and Hammer, 1993). Unfortunately, this theoretically 
predicted axial ratio distribution has not been seen in the recent arc survey with the subsample of EMSS 
clusters. Among the 14 arcs found in the 41 EMSS clusters, 9 arc giant arcs (Hammer et al., 1995; Luppino 
et al., 1995). It is believed that the asymmetrical matter distribution of the arc clusters may account for the 
discrepancy (Bartelmann and Weiss, 1994; Bartelmann, Steinmetz and Weiss, 1995). 

Modeling the known giant luminous arcs, even arclets, associated with clusters of galaxies turns to be 
very successful. A simple elliptical or bimodal potential for the arc cluster can provide the major properties 
compatible with the observed arcs. Indeed, any type of elongated configurations including the straight arcs 
and the radial arcs has been well reproduced based on elongated gravitational potentials tracing the luminous 
matter distributions (Pello et al., 1991; Mellier, Fort and Kneib 1993). An example of modeling the arcs in 
the well-known Abell cluster 370 is shown in Figure 30. This arc-cluster system has been extensively studied 
by a number of authors (e.g. Hammer and Rigaut, 1989; Grossman and Narayan, 1989), in particularly by 
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Figure 30: Modeling the giant luminous arc and arclets in Abell 370. Upper left: The direct CCD image of 
the core of Abell 370; Upper right: Image construction. Solid lines and dashed lines are, respectively, the 
caustics and critical lines (the inner one is for Zg = 0.725 and the outer one, Zg = 0.895). Dot-dashed lines 
are the core radius, ellipticity and orientation of the two matter profiles: Low left: Contours of the surface 
density; Low right: Reconstruction of the source positions (see the upper right panel for the corresponding 
part), (from Kneib et al., 1993) 



the Toulouse Group (e.g. Kneib et al., 1993; Soucail and Mellier, 1993). 



4.3 Cluster matter distributions from arcs 

Arclike images are robust matter estimators of clusters of galaxies, which is actually the most important 
issue of studying arclike images today. Prom the general expression of the lensing equation [eq.(12)] for a 
spherical matter distribution [eq.(8)], we can write out the projected mass along the line of sight within the 
arc position 9 ( in arcseconds) to be 

mg{e) = 7.37 X 10^\e - l3)e^^^B^ Mq /igoS (100) 

where /3 is the alignment parameter of background galaxy in arcseconds. This parameter, however, is 
unmeasurable in practice. Nevertheless, if we further assume a SIS model for the matter distribution of 
cluster of galaxies and a uniform circular disk for background source, the max:imum width (W) of the arclike 
image will be the same as the size of the background galaxy and /3 will satisfies the following approximate 
geometrical relation 

2{L/W)sm{L/2eDd)' 

Therefore, the gravitational mass contained in the central core (< ^) of an arc cluster is obtained without any 

assumptions about the dynamical state of the system. For most of the arc-cluster systems, the gravitational 
mass is typically ~ lO^'^M© in the core of a galaxy cluster. The significance of this method is that it provides 
an independent way of calculating the masses of clusters of galaxies, which can be compared directly with the 
masses estimated from the dynamical analysis based on the virial theorem. Recall that the latter presumes 
an hydrostatic equilibrium of both the hot intracluster gas and the galaxies with the binding cluster potential 
(Cowie, Henriksen and Mushotzky, 1987). This hypothesis, unfortunately, has not been verified using other 
astrophysical means. 

Assuming that both the X-ray gas and the galaxies are in hydrostatic equilibrium with the binding 
gravitational potential of a spherical galaxy cluster, one can obtain the virial mass of the cluster within 
radius of r through 

kTr /dlnn dlnT' 



Mv{r) = --^ -TT— + -TT— , (102) 

where Hp is the mean particle weight in unit of the proton mass m,p, n and T are the gas density and 
temperature, respectively, which can be found by inverting the observed X-ray surface brightness profile, 
namely, the /? model 

5(^) = 5o[l + (W']'^'"''^ (103) 
with a core radius of 0c (or Tc in linear size) . The isothermal gas distribution has been found to be consistent 
with the X-ray observations of galaxy clusters, leading to dlnT/dlnr = 0. So, the projected virial mass 
my{0) from eq.(102) within a radius of on the cluster plane is (Wu, 1994c) 

m.{6) = 1.14 X 10"^m(^) M^; (104) 
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Table 6: Four arc-cluster systems and their masses 



rh{9) = 



l + Rl 



Rl 




91 (1 + ,t2) 



-dx. 



(105) 



Here i?o = R/tc, and R is the physical size of the cluster. The numerical computations show that m{9) 
remains nearly unchanged for R ranging from 3 Mpc to 100 Mpc. 

A comparison of \\iv gravitational masses estimated from arclikc images [cq.(lOO)] with the masses derived 
from hydrostatic equilibrium [eq.(104)] is shown in Table 6 for four clusters of galaxies in which both arclike 
images are detected and the X-ray data are available (Henry et al., 1982; Gioia and Luppino, 1994): Abell 
370 (A5), MS 1006.0+1202 (arc 4), MS 1008.1-1224 (arc 2) and MS 1910.5+6736. It appears that there 
exists a significant difference between the virial masses and the gravitational masses from arclike images in 
all the four systems. 

Three groups have independently announced similar results for a total of 8 arc-cluster systems (Wu, 
1994c; Fahlman et al., 1994; Miralda-Escude and Babul, 1995), showing that the virial equilibrium has 
underestimated the total gravitational masses of clusters of galaxies by a factor of at least 2.5 up to the arc 
positions. There are three main possibilities that may account for the mass discrepancy: (l)the hot gas in 
clusters of galaxies may be meanwhile supported by a non-thermal pressure such as magnetic field (Loeb 
and Mao, 1994; En/?lin et al., 1995); (2)cluster matter distributions may be highly prolate with the long axis 
along the line of sight (Miralda-Escude and Babul, 1995); And (3)clusters of galaxies cannot be considered 
to be the well relaxed virialized systems (Wu, 1994c). Nonetheless, the third possibility, if true, may offer 
an important clue to resolving the "baryon catastrophe" on scale of clusters of galaxies (White et al., 1993). 

The cosmic baryon fraction is the ratio of baryonic matter Mf, (X-ray gas + galaxies) to the total mass 
M (baryon + non-baryon) of clusters of galaxies: = Mi,/M, provided that clusters of galaxies are repre- 
sentative of the matter distribution of the Universe. The baryon catastrophe arises because the virial mass 
My derived from eq.(102) is used as the measurement of the total mass M of cluster of galaxies, which has 
led to rib being 3 ~ 10 times larger than the prediction of the Big Bang Nucleosysthesis and the standard 
inflation cosmological model. Now, replacing the virial mass my by the gravitational mass rUg estimated 
from arclike images would reduce the baryon fraction ^5 by a corresponding factor of at least 2.5 over the 
region of the typical core radius of cluster. As a consequence, the "fib discrepancy problem" might vanish. 



4.4 Cluster matter distributions from statistical lensing 

Another useful constraint on the matter distribution of galaxy clusters is provided by the study of statistical 
lensing, which attempts to statistically investigate the possible form of matter distributions of galaxy clusters 
as a whole using the properties of giant arcs and/or arclets, such as the total number, the width and the 
axial ratio of arcs (Hammer, 1991; Wu and Hammer, 1993; Miralda-Escude, 1993a, b; Grossman and Saha, 
1994). This procedure is actually a convolution of the magnification probability of galaxy clusters with the 
distribution of background galaxies. 

The differential probability that a source at Zs is magnified by a factor within of /i due to a galaxy 
cluster at Zd is proportional to the lensing cross-section 27r£)^/3d/3 [eq.(ll)]. To find the total magnification 
probability P{zs, n) even for a spherical matter distribution of cluster of galaxies, one still needs to know the 
cluster velocity dispersion distribution, the length scale distribution and the spatial distribution. Unfortu- 
nately, all these distributions have not been very well determined from observations. The X-ray observations 
of clusters of galaxies from the Einstein Observatory and EXOSAT (Edge et al., 1990) as well as the EMSS 
(Gioia et al., 1990; Henry et al., 1992) provide a Schechter luminosity function at low rcdshift {zd < 0.2) and 
a power-law luminosity function <l)[Li4)dLx = KL^^dL^ [^44 = {Lx/IQ'^^ ergs s""'^)] at intermediate redshift 
ranging from 0.14 to 0.6, showing a significant evolution of X-ray luminous clusters with cosmic epoch, where 
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Figure 31: Magnification probability of gravitational Icnsing by clusters of galaxies for four density profiles. 
The variable rc distribution of Jones and Forman (1984) is adopted for ISC and KING models and a constant 
Te of 6h^o Mpc is used for the r^/* law. A no-evolution scenario is assumed for the number distribution of 
galaxy clusters. Furthermore, the extension of the background source is neglected. 



K takes different values at different redshift shells and has units of Mpc~^[I/44]'^~^, and v is the power-law 
index. The X-ray luminosity distribution can be converted into the velocity dispersion distribution through 
the correlation between X-ray luminosity and velocity dispersion (Quintana and Melnick, 1982; Wu 
and Hammer, 1993): 

i, = 1032-7V^-94ergss-i. (106) 

For the distribution of length scale (core radius rc, effective radius re, etc.), optical observations exhibit a 

constant core radius of rc = 0.25h'^Q Mpc for the King model (Bahcall, 1975;1977) and a constant effective 
radius of re = 6h^o Mpc for the r^^^ law (Bears and Tonry, 1986), while X-ray observations provide a 
distribution of core radius significantly different from the constant core (Jones and Forman, 1984): 

Pe(logre)dlogre = ^ cxp [- . 5 (log re + 0.89^)/0.32^] d log re, (107) 
0.32v27r 



which leads to an average core radius of O.llh'^Q Mpc. Both the constant and variable scale length dis- 
tributions will be adopted in the computation of P(zs,/x). As for the number density of galaxy clusters, 
both no-evolution and evolution models from the X-ray data will be used and the respective results will be 
compared. 

The total magnification probability P{zs, jj) for a source at = 2 by intervening clusters of galaxies with 
different matter distributions is shown in Figure 31. Significant differences in P{za, jj) due to different models 
of mass density are clearly seen. This property opens then a possibility to test the matter distribution of 
clusters of galaxies using statistical lensing. Convolving P{zs, ji) with the distribution of background galaxies 
would yield the number of lensed galaxies with magnification greater than fi. However, n measures only 
the magnification of the apparent luminosity rather than the geometrical features of the lensed source. The 
relation between the magnification ji and the axial ratio L/W oi the elongated image of a background circular 
galaxy with radius of Rq by a spherical lens can be approximately obtained to be (Wu and Hammer, 1993) 



L 



l+2)o4?^-27ri5oS(^o) 

^0 



X ^sm (108) 



in which we have presumed that W/Oq <^ 1, i.e., the arc width (W) is much smaller than the arc distance (^o) 

from the center of its associated cluster (Apparently, many detected arcs can meet this condition) , where /3o 
and 6*0 are the alignment parameter and the corresponding image separation for the center of the background 
galaxy, and Dq = {4:G /c^){DciDcis / Dg) . In particular, as giant arcs trace the critical line or the Einstein ring 
9e, do can be approximately replaced by 9e- Furthermore, one can expand sin~"'^(i?o//3o) in (-Ro/Z^o) if the 
complete ring images are excluded in the statistics due to their rareness. Then, eq.(108) becomes 

A=4Ml-Kf(l-if + ...), (109, 

where K = 'k6%D'^Tj{9 e) / ^^{O e) ■ Numerical computations using the three spherical models in Table 1 show 
that the term (l/6)(i?o//3o)^ becomes important only if the radius of the circular galaxy is larger than 8 kpc 
and jjL > 30, while most of the known giant arcs have ^ 10. Therefore, neglecting the source extension 
cannot cause too serious problems in the evaluation of P{zs, ii). Thus, 

.= ^^. (110) 
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Table 7: The Expected Number of Giant Luminous Arcs over the Whole Sky 
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This expression relates magnification with axial ratio of the distorted image and is also very useful to estimate 
the image magnification from its geometrical configuration. For example, for the SIS a simple calculation 
shows K — 1/2, which then leads to /x = L/W. That is to say, the magnification remains the same as 
the axial ratio for the image produced by SIS and thereby, the magnification probability P{zs,^i) is the 
probability that a source is elongated by a factor of greater than ji = L/W. For other spherical mass density 
models one can find the probability P{zs,^, L/W) from eq.(llO) that a source at Zs is magnified by a factor 
of fi and elongated by a factor oi L/W. 

The expected number of giant arcs is finally obtained by convolving the probability P{zs,fi, L/W) with 
the Schechter luminosity function of background galaxies which is taken from the survey of Broadhurst, Ellis 
and Shanks (1989). Table 7 gives the result of Wu and Hammer (1993) for the totally expected number of 
giant luminous arcs with B < 22.5 and L/W > 10 within Zs < 1.25 over the whole sky. Four mass density 
models arc used for the matter distribution of clusters of galaxies which are restricted within the range of 
0.15 < Zd < 0.6. Both the no-evolution [the Edge et al. (1990) local X-ray Schechter luminosity function] 
and the evolution model [the Henry et al. (1992) X-ray power-law luminosity function] of galaxy clusters are 
used. As a consequence of the significant difference in P{zs,fj,) arising from the various matter distributions 
(see Figure 31), the resulting number of giant luminous arcs shows remarkable differences. The strong 
evolutionary scenario of the X-ray luminous clusters yields fewer clusters than does the no-evolution model, 
giving rise to fewer giant luminous arcs predicted from the evolution model of clusters. Apparently, very rich 
clusters of <t„ > 1300 km/s have a relatively higher frequency of producing giant arcs than the rich clusters 
of av > 800 km/s. A further computation using the clusters of cr„ > 1300 km/s within the redshift interval 
0.15 < Zrf < 0.4 gives a frequency of detecting giant luminous arcs to be as high as 30-40%. Nevertheless, 
the model ISC in Table 7 has underestimated the number of giant luminous arcs by a factor of at least 2 
as compared with the known giant luminous arcs unless the core radii of clusters are reduced significantly, 
indicative of a smaller core for the total matter distribution than the X-ray gas. In the extreme case of the 
zero-core radius, SIS predicts a number of giant luminous arcs that marginally reconciles with the observed 
arcs in the strong evolution scenario of clusters. It is then concluded that the well-known form of ISC, which 
is found to be the best-fit to the X-ray luminosity distribution in clusters of galaxies, cannot account for the 
giant luminous arcs if dark matter has the same core radius as the observed one in X-ray. This implies a 
more compact dark matter distribution in the central core of the galaxy cluster. 

The statistical study of arc widths has also reached a similar conclusion (Grossman and Saha, 1994) 
that cluster mass density profiles must have core radii at least as small as 0c/Oe < 0.1 in ISC. Interestingly, 
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the dynamical analysis of a set of 12 clusters of galaxies using the X-ray data from the Einstein Imaging 
Proportional Counter has found that dark matter appears more peaked in the cluster centers than the X-ray 
gas and has core radii of only 50—100 kpc (Gerbal et al., 1992; Durret et al., 1994). Although this dynamical 
method has given a result that agrees with the lensing estimate, it is not so clear if the dynamical masses 
derived from the assumption of hydrostatic equilibrium for galaxy clusters are reliable values in view of the 
argument in section 4.3. 

It must be pointed out that the above conclusions drawn from statistical lensing by clusters of galaxies 
should be taken to be very preliminary. Wu (1993d) has discussed four parameters in the statistical lensing 
which may cause some large uncertainties in the predictions of arc properties due to their limited knowledge 
from today's observations: length scale, velocity dispersion, extended source and evolutionary effect. (l)Core 
or effective radii: Both optical and X-ray data are not sufficient to constrain the distribution of core/effective 
radii of galaxy clusters. The optical core radii claimed by Bahcall (1975, 1977) are two times smaller than 
those fomid by Dressier (1978), whilst the distribution of X-ray core radii eq.(107) shows somewhat larger 
scatters. Alternatively, the constant effective radius of Qh^Q Mpc adopted for the r^^^ law is apparently not 
a good value for different rich clusters of galaxies. So, the arc predictions based on the poor data of length 
scale of clusters of galaxies may have large imcertainties. (2)Velocity dispersion: Magnification probability 
for lensing is very sensitive to this parameter. But the average line of sight velocity dispersion along the 
radius of a galaxy cluster is a variable rather than a constant except for SIS (Wu, 1993b; Kochanek, 1993c). 
The above statistical lensing takes into account only the asymptotic value of cr^, , either close to the center or 
at large radius, which may result in a significant variation in lensing probability. (3)Extended sources: The 
above point source approximation for background galaxies would overestimate the magnification probability 
P{zg,fi). Recall that a maximum magnification iimax instead of infinity is reached for an extended source 
(see section 1.3). Moreover, the consideration of source extensions might solve the puzzle of too many arclets 
relative to the giant arcs predicted from statistical lensing (Bergmann and Petrosian, 1993). (4)Evolution- 
ary effect: The strong evolution of galaxy clusters with cosmic epoch would lead to a number decrease of 
clusters with redshift, providing much fewer lenses at high redshift for background galaxies. Conversely, the 
no-evolution model is able to produce a relatively larger number of giant luminous arcs in the clusters of 
galaxies at high redshift. The observed giant luminous arcs then open a possibility of testing the evolution 
model of distant clusters of galaxies (Wu, 1993a; Bartelmann and Weiss, 1994). Nevertheless, the poorly es- 
tablished evolutionary model of clusters of galaxies may affect the present computation of lensing probability. 



4.5 Cluster matter distribution from weak lensing 

The parametric likelihood method of determining the cluster matter distribution employed in sections 4.3 
and 4.4 is suitable for modeling the giant arcs/arclets, for constraining the length and shape of global cluster 
potential wells and for investigating the relative redshift distribution of the background objects. The method 
assumes a priori a density profile of the cluster and determines the most likely values of the model parameters 
using the properties of the observed arcs/arclets. It appears that the current data of the distorted images 
are still insufficient to discriminate between the various models. It is, therefore, most desirable that one can 
directly derive the cluster matter distribution with a parameter-free (or non-parametric) method. 

Tyson, Valdes and Wenk (1990) firstly detected the coherent distortion (weak lensing) of faint blue 
galaxies behind two clusters of galaxies (A1689 and CL 1409+52). They used the alignment statistics to 
extract the lensing signal which is characterized by an ellipticity (Valdes, Tyson and Jarvis, 1983) 

(111) 



a2 + 62 ' 

where a and b are the semimajor and semiminor axes of the principal axis transformed moments. This 
yields a positive value of the net alignment for a population of background galaxies around a foreground 
cluster, in comparison with the zero result for a random population of galaxies. Kaiser and Squires (1993) 
extended this idea to the outer parts of clusters and to statistical cluster samples, and developed a powerful 
technique to reconstruct the cluster surface mass density S(0) through the measured distortion e, which 
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is totally independent of the presumed density profile for clusters. This non-parametric method has been 
successfully used in the mapping of two-dimensional mass distributions in several clusters and is now being 
further developed for its wide applications in both weak and strong lensing situations (Seitz and Schneider, 
1994; Kaiser, 1995; Bartelmann, 1995: Bartelmann and Narayan, 1995; etc.). 

The image shapes are characterized by the quadrupole moments (Valdes, Tyson and Jarvis, 1983) 



/ (fee,eji{e) 
JdmW 



Qij = ' VZua. (112) 



jdm'id) ^ "^^ 

where I{0) and I' (9) are the intrinsic and observed surface brightness distributions of a background source 
[see also eq.(17)], and the angles are measured from the centroid of the image. Assuming that the source 
is relatively small so that the magnification matrix A is constant over the images, we have the following 
transformation between the quadrupole matrix of source and image according to eq.(17) 

Q,j = AikAjeQ'kt, (114) 



or 



Q = AQ'A, (115) 
where A is defined hy A= dfi/dO, and in terms of eqs.(18) and (19) 
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(116) 



By analogy with eq.(lll), we define the intrinsic ellipticity parameter e of the background source and the 
observed ellipticity parameter e' of the image as 

e' = {e[,e',) = (%^^, J^''' ], (118) 



k Q'li "I" Q22 Q'li "I" Q22 J 

which can also be denoted in their complex forms: e = ei + ie2 and e' = e'i+ ie2. 

In the limit of weak lensing, one can use the linear approximation. Under the transformation eq.(115), 
eqs.(117) and (118) reduce to 

e[ = ei + (V'li - V'22)(l - ef) - 2V'i2eie2; (119) 
62 = 62 - (V'li - V'22)eie2 + 2^/;i2(l - el). (120) 

The mean intrinsic ellipticity of an ensemble of background galaxies should be zero: (ei) = (62) = (6162) = 0, 

while the factors 1 — {ef) and 1 — (e|) are close to unity in practice. As a result, the expectation of the image 
ellipticity is simply 

(e'l) = Vii - V'22; (121) 

{e'2) = 2V;i2, (122) 

or 

(e^) = 27i, (123) 

i.e., the image distortion is uniquely determined by the tidal field 7 of the foreground clusters. However, the 
shear term 7 is related to the matter term k through [eqs.(14). (16) and (19)] 

ji0) = ^Jd^e'x{e-d')K{e'), (i24) 
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Table 8: Quasar-cluster associations: observations and models 



where the kernel x(^) is a complex function (Schneider and Seitz, 1995) 

xiO) = (125) 

Finally, the inversion of eq.(124) yields the cluster surface mass density S(0), which is expressed by the 
observed ellipticity parameter of the images, 

^(^) = V / ^'^'^e[x*(0 - e')e'{e% (126) 

where x* is the complex conjugation of x and Re takes the real part of the complex variable. 

Since the first detection of gravitational weak shear out to a radius of 3.0 h'^Q Mpc in CL 0024+1654 
(Bonnet, Mellier and Fort, 1994), several measurements of the weak distortion of background galaxies have 
been made and the applications of the lensing inversion technique to the observed data have turned to be very 
successful in the reconstruction of the mass distribution of clusters (Fahlman et al., 1994; Small et al. 1995; 
Small and Dickinson, 1995; Tyson and Fischer, 1995; Squires et al., 1995; Kneib et al., 1995; etc.). Although 
there are some disagreements about the resulting gravitating masses of clusters, for instance, Fahlman et 
al. (1994) derived a gravitational mass of 2.5 - 3 times larger than its virial mass in cluster MS 1224 while 
Squires et al. (1995) found accordance between the two masses in Abell 2218, it appears very promising 
that one can precisely determine the cluster mass from lensing method in the near future by improving the 
inversion technique (Bartelmann, 1995; Schneider, 1995; Schneider and Seitz, 1995; Seitz and Schneider, 
1995a; Kaiser, 1995). 



4.6 Quasar-cluster associations 

The presence of giant arcs and arclets associated with gravitational potentials of clusters of galaxies indicates 
that clusters of galaxies are very efficient lenses. Actually, it was noticed at the same time when giant 
luminous arcs were discovered that some of the high-redshift 3CR galaxies are gravitationally magnified by 
the low-rcdshift clusters of galaxies lying in their lines of sight (Hammer, Nottale and Le Fcvre, 1986; Le 
Fevre, Hammer and Jones, 1988; Le Fevre et al., 1988; Hammer and Le Fevre, 1990). Motivated by these 
observations, Wu and Hammer (1993) even explored the possibility of whether there are radio arcs behind 
galaxy clusters. 

What would happen to background quasars if foreground clusters of galaxies act as lenses ? Four recent 
measurements have answered this question by discovering a significant quasar overdensity behind foreground 
clusters using different quasar and cluster samples: (l)thc Large Bright Quasar Survey and Zwicky clusters 
(Rodrigues- Williams and Hogan, 1994), (2)the 1 and 2 Jy Radio Source Surveys and Abell clusters (Wu and 
Han, 1995), (3)the variability selected quasars and clusters (Rodrigues- William and Hawkins, 1995). (4)the 
1 Jy Radio Source Catalog and Zwicky clusters (Seitz and Schneider, 1995b). Table 8 summarizes these 
searches and their resulted quasar overdensity density q, in which only the most significant q measured at 
a fixed position 6 and a limiting magnitude (or flux) is given. Figure 32 shows the variations of q for Abell 
clusters at Zd < 0.2 aroimd 2 Jy radio quasars at Zg > 0.5. 

In a similar way to the study of quasar-galaxy associations, we can evaluate the enhancement factor q for 
the quasar-cluster associations. Given a magnification /x, a general expression is available for both optically- 
and radio-selected quasars 

^ iV(<m + 2.51ogM) 1 ^ ^(> S/ii) 1 
* N{<m) II N{>S) n ^ ' 
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Figure 32: Variations of enhancement factor of Abell clusters around 2 Jy radio quasars. The data have been 
normaUzed at 3° and the standard deviation (la) is plotted in each measurement. The fit to the prediction 
by a singular isothermal sphere as lensing object is shown using a critical radius of 0.2°. 



where m and S denote the limiting magnitude and the flux threshold, respectively. The surface number 
density of optically-selected quasars has been given by BSP [see eq.(94)]. For the radio source counts 
N{> S), a least-square fit of a power-law to the observations at 5 GHz (Langston et al., 1990; Fomalont et 
al., 1991) yields 

^(>'^)- I 2.10 X 105 5-110, 5<10mJy, ^^^^^ 

where S is the units of mJy. Note that radio source catalogs are composed not only of galaxies but also 
quasars. The fraction of quasars in radio source surveys varies with flux threshold. Therefore, the employ- 
ment of N{> S) in the study of quasar-cluster associations provides only an estimate of q. 

Now we work with the lensing models of clusters of galaxies and/or their associated matter inhomo- 
geneities and test whether one can explain the reported quasar overdensity behind clusters on scale of ~ 10 
arcminutes in terms of gravitational lensing. We take an average enhancement (q) [eq.(93)] over the search 
range of 9 around clusters instead of the local enhancement q at d. We further assume a flat cosmological 
model n — 1 and adopt Hq = 50 km/s/Mpc. 

Conventionally, clusters of galaxies should be considered to be the lensing objects for the reported quasar- 
cluster associations. Utilizing SIS as the mass model and its magnification of eq.(27), we can estimate the 
cluster velocity dispersicm (T,; that is required to produce the observed (q). Figure 32 plots such a fit to 
the observed data. Surprisingly, the best fitted Einstein radius in this example is ^ 0.2°, corresponding to 
a velocity dispersion of tXc ~ 5000 km/s if the typical redshifts of Abell clusters and of radio sources are 
taken to be 0.1 and 1, respectively. The similar results are found for the rest three measurements (Table 
8). Apparently, the masses (~ a^) that are needed to produce the four measured enhancement factors are 
substantially larger than the realistic value for clusters. 

It has been known for some years that the weak lensing by large-scale matter inhomogeneities may con- 
tribute a significant effect on the background quasars. It may be the cause for the quasar-galaxy associations 
observed on the similar scale (~ 10') (Fugmann, 1988;1990; Bartelmann and Schneider, 1993a,b;1994). We 
can now work out how largc^ an additional mass surface density from the large-scale matter chimps that 
clusters of galaxies trace is need to produce the quasar-cluster associations. To do this, we add a uniform 
mass sheet S to clusters of galaxies. It turns out that the Einstein radius 9e and the image separation are 
increased by a factor of (1 — E/Ec)-i and the lensing magnification becomes 



1 

\i-eE/e\(i-E/Eey 



/^= H n_,nu. ^/^^2 > (129) 



where Sc is the critical surface mass density that any lens must exceed in order to produce multiple images 
by itself [eq.(15)]. Quantitatively, the minimum Sc for a source at = 2 is 0.41 g/cm^. In Table 8 we 
give the required surface mass density S for each of the measurements. Note that S deduced from the radio 
selected quasars is a factor of ~ 2 larger than the one from the optically selected samples. This is due to 
the contamination of radio galaxies in the radio source catalog. The result from Seitz and Schneider (1995) 
illustrates very well this effect. 

We can estimate the matter contribution from all the galaxy clusters that follow the cluster spatial two- 
point correlation function ^{r/rcc}^^'^ ■ where the correlation amplitude is = 40 Mpc (Postman, Huchra 
and Geller, 1992). If we assume that each cluster has a SIS mass density profile and a gravitational radius 
Rc and the mean number density of clusters is constant, the surface mass density of clusters enclosed within 
6 around a given cluster at Zd is (Wu and Fang, 1995) 

S(^) = 4noMcrcc(l + zaf F{e, Vcc, Rc), (130) 
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where no is the present cluster number density, Mc = 2a'^Rc/G is the total cluster mass and F is a function 
given by $, (r) and SIS model. Numerical computations show that « 2 ~ 3 for i?c = 3-5 Mpc over the range 
of ^ = l'-80'. Let flc denote the fraction of the total cluster matter in the mass density of the Universe, we 
have 

I] = O.Ome(^^^)(f) g/cm^ (131) 

Unfortunately, this surface mass density provided by all galaxy clusters following ^(r) is an order of magnitude 
lower than that required for the quasar-cluster associations even if Qc = 1. 

The mass surface density from large-scale structures of the Universe can be estimated through 



/ [Pir) 



po]dr ~ 1.45 X 10--'S { j^^^ ] g/cm^ (132) 



in which S is the mean density contrast over scale of R. However, the evaluation of S is sharply constrained 
by the measurements of temperature anisotropy AT/T of the cosmic background radiation on various scales. 
Numerical computation indicates that it is impossible to attribute the large mass surface density derived 
from the quasar-cluster associations to any matter clumps on scale of i? > 20 Mpc in the Universe if 
AT/T = 1 ~ 5 X 10-5 (Wu and Fang, 1995). 

So, if the reported associations between background bright quasars and foreground clusters of galaxies are 
not due to statistical variations arising from the quasar /cluster selections and patchy Galactic obscuration, 
we need to consider the following possibilities: (l)There may exist a large amount of unseen matter between 
clusters of galaxies on scale of ^ 10 Mpc, because the above calculations did not include the unbound clus- 
ter matter. This can be tested using the N-body simulations, as was made by Bartelmann and Schneider 
(1993b) for quasar-galaxy associations. (2)The working hypothesis may be wrong, i.e., the observed back- 
groimd quasar coimts may deviate from their intrinsic ones. On the scale of galaxies, Schneider (1992) has 
demonstrated that dropping the unaffected background hypothesis does not significantly improve the situa- 
tion in quasar-galaxy associations. Whether the cluster matter or large-scale structures would contribute a 
non-negligible effect on the quasar number counts needs to be further investigated. 



FINAL REMARKS 

The past few years have been exciting times for lensing people. In particular, the microlcnsing experiments 
have detected a few ten events associated with the compact objects of the Galactic halo/disk and/or the 
LMC halo/disk, indicative of the success of using gravitational lensing effect for the searches of dark matter, 
which will have a strong impact on various aspects of cosmology study today. The rapidly increasing 
number of new lens systems (multiple quasars/galaxies, radio rings/galaxies, arcs/clusters of galaxies, etc.) 
has made it possible to study the matter distribution of the Universe statistically and to determine the 
cosmological parameters {Hq, and Ao). This is of particularly significance since lensing provides not only 
an independent means to evaluate these important issues in cosmology but also a test for the validity of 
other astronomical/physical methods. The determination of the Hubble constant Hq from the time delay of 
double quasars and mapping matter distribution in clusters of galaxies with luminous arcs and arclets are 
the two excellent examples of the lensing applications in cosmology. 

Indeed, the study of gravitational lensing has been developed so rapidly in both theory and observation, 
and it is not practical to summarize every subject in this review destined for the readers who are not 
the experts in lensing. In the present article it is even impossible to include some new discoveries, new 
observations and new theories which appeared during the writing of the article. The interested readers 
are recommended to refer to the recent reviews on the lensing applications in cosmology (Blandford and 
Narayan, 1992; Schneider, 1996), on the lensing observations (Refsdal and Surdej, 1994), on the arc(let)s in 
clusters of galaxies (Fort and Mellier, 1994), in particular, the excellent monograph of gravitational lensing 
by Schneider, Ehlers and Falco (1992). 



39 



Two quotes can be used as the final remarks of this review on gravitational lensing: 

"An astronomer can use beams of photons to probe a condensation of dark matter in much the same 
way that a nuclear physicist used beams of electrons to study the structure of an atomic nucleus. " (Bland- 
ford and Kochanek, 1987) 

"A galactic gravitational lens can be used as the ultimate astronomical telescope. " (McBreen and Met- 
calfe, 1987) 
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Table 1 Gravitational Lensing Models 



model 


ISC 


KING 


r^/* law 


surface density 


So 


So 


Soexp[-7.6696'o^*] 


length scale 


rc 




re 


central density (po) 
or total mass (Mo) & 
velocity dispersion (ay) 


PO — 27rGr§ 


PO — i-KGrl 


Mo = 9.0^ 


lensing equation 






l3o = eo-D^ (*) 


D 




C' rcDs 


364 


critical D 


2 


1 


3.370 X 10-3 



* mo(^o) = 1 - exp (-7.6696*0^^) ^ ^ (jMdOl^'^^ 
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Table 3 Gravitationally-Lensed Multiple Quasars* 



name 


image No. 


lens 




Zs 




status 


discovers 


year 


0957+561 


2 


G+G 


0.36,0.5 


1.41 


6 .1 


coniirmed 


Walsh et al. 


1979 


1115+080 


> 4 


G 


0.29 


1.72 


o// o 

2 .3 


coniirmed 


Weymann et al. 


1980 


2345+007 


2 


G(?) 


1.49(r) 


2.15 


7 .3 


possible 


Weedman et al. 


1982 


iDo4+ZD/ 


o 
A 


<^(-) 




1.90 


Q" C 
O .0 


possible 


Djorgovski & Spinrad 


19o4 


zOlo+lli 


o 

3 


G 


1.01(?) 


O OT 

3.27 


3 .8 


connrmed 


Lawrence et al. 


1 no A 

1984 


2237+0305 


4 


G 


0.04 


1.69 


1 .8 


connrmed. 


TT 1 J. 1 

Huchra et al. 


1985 


(J142-1UU 


2 


G 


0.49 


O TO 


o// o 
2 .2 


coniirmed 


Surdej et al. 


1987 


1413+117 


4 


G(?) 


1.4(?) 


2.00 


1.1 


coniirmed 


Magain et al. 


1 r\oo 

1988 


1 1 on 1 m n 
iizU+Oiy 


o 


G(OG(0 


0-6(0 


i.40 


.0 


possible 


Meylan & Djorgovski 


1 OQO 
LVOV 


1429-008 


2 


(?) 


1.6(?) 


2.08 


5.1 


possible 


Hewett et al. 


1989 


nn CO ni 1 c 
0952-0115 


2 


v) 


(-0 


4.5 


U .9 


possible 


McManon et al. 


1 nno 
1992 


1208+1011 


2 


(?) 


(?) 


3.80 


0".47 


possible 


Magain et al. 


1992 


1422+231 


4 


G(?) 


0.64(?) 


3.62 


1".3 


coniirmed 


Patnaik et al. 


1992 


1009-025 


2 


G(?) 


1.62(?) 


2.74 


1".55 


possible 


Surdej et al. 


1993 


1104-1805 


2 


(?) 


1.66(?) 


2.30 


3" 


possible 


Wisotzki et al. 


1993 


0240-343 


2 


G(?) 


0.34(?) 


1.4 


6".l 


possible 


Tinney 


1995 


J03.13 


2 


G(?) 


1.09(?) 


2.55 


0".84 


possible 


Claeskens et al. 


1995 



* Radio sources are not included. G=galaxy; C=cluster; A^=separation 
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Table 5 Arclike Images 



arc cluster Z4 cr^(km/s) Z/x,44(ergs/s) arc redshift(2;s) discovers year 



AbcU 222 


0.213 


570 


3.7 




Small et al. 


1991 


Abell 370 


0.374 


1364 


9.7 


0.725,1.3? 


Soucail et al. 


1987a 


Abell 963 


0.206 




9.1 


0.771 


Lavery & Henry 


1989 


AbcU 1689 


0.196 


1989 


17. 




Tyson et al. 


1990 


Abell 1942 


0.224 








Small et al. 


1991 


Abell 2104 


0.155 




8.0 




Pierre et al. 


1994 


AbcU 2163 


0.203 






0.728,0.742 


Soucail ct al. 


1994 


Abell 2218 


0.176 




6.5 


0.702,1.034 


Fello-Descayre et al. 


1988 


Abell 2219 


0.225 




18. 


~ 1 


Small et al. 


1995 


AbcU 2280 


0.326 




5.1 




Gioia ct al. 


1995 


Abell 2390 


0.231 






0.913 


Fello-Descayre et al. 


1991 


Abell 2397 


0.212 








Small et al. 


1991 


Abell S295 


0.301 


900 






Edge et al. 


1994 


CI 0024+1654 


0.391 


1300 


2.7 


1.39? 


Koo 


1987 


CI 0302+1658 


0.426 




5.0 




Mathez et al. 


1992 


CI 0500-24 


0.316 


1375 




0.91? 


Giraud 


1988 


CI 1409+52 


0.46 


3000 


9.2 




Tyson et al. 


1990 


CI 2236-04 


0.56 






1.116 


Melnick et al. 


1993 


CI 2244-02 


0.336 




1.5 


2.237 


Soucail et al. 


1987a 


MS 0440+0204 


0.190 




4.0 




Luppino et al. 


1993 


MS 0451-0305 


0.55 




20. 




Lc Fcvrc ct al. 


1994 


MS 1006+1202 


0.221 




4.8 




Le Fevre et al. 


1994 


MS 1008-1224 


0.301 




4.5 




Le Fevre et al. 


1994 


MS 1455+2232 


0.259 




16. 




Lc Fcvrc ct al. 


1994 


MS 1621+2640 


0.426 




4.5 




Le Fevre et al. 


1994 


MS 1910+6736 


0.246 




4.4 




Le Fevre et al. 


1994 


MS 2053-0449 


0.583 




5.8 




Le Fevre et al. 


1994 


MS 2137-2353 


0.313 




16. 




Fort et al. 


1992 


MS 2318-2328 


0.187 




6.8 




Le Fevre et al. 


1994 


AC 114 


0.31 


1649 


4.0 


0.639 


Small et al. 


1991 


0956+561* 


0.36,0.5 








Bernstein et al. 


1993 


GHO 2154+0508 


0.32 






0.721 


Lavery et al. 


1993 


PKS0745-191 


0.1028 






0.433 


Allen et al. 


1995 


RXJ1347.5-1145 


0.451 




62. 




Schindler et al. 


1995 


ACQ 3408 


0.042 






0.073 


Campusano & Hardy 


1995 



* The recent observation (Dahle, Maddox and Lilje, 1994) suggests that this "arc system" is the result of chance 
aligments of three and two different objects, and not gravitationally lensed arcs. 
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Table 6 Four arc-cluster systems and their masses 



cluster 


arc 


mass 


name 


Zd 




en 


L{") 


L/W 


Zs 


m^{e) (Mq) 


mg{e) (Mo) 




Abell 370 


0.374 


9.7 


56 


9 


18 


1.3(?) 


2.26 • 10" 


8.20 • 10" 


3.63 


MS1006.0+1202 


0.221 


4.819 


62 


4.9 


7.0 


0.6 
2.0 


1.36- 10" 
1.36- 10" 


6.91 • 10" 
4.87- 10" 


5.08 
3.58 


MS1008.1-1224 


0.301 


4.493 


51 


4.0 


6.5 


0.6 
2.0 


1.33- 10" 
1.33-10" 


7.47 • 10" 
4.34 • 10" 


5.60 
3.25 


MS1910.5+6737 


0.246 


4.386 


67 


6.1 


10.5 


0.6 
2.0 


1.63- 10" 

1.63- 10" 


9.95 • 10" 

6.64- 10" 


6.12 

4.08 
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Table 8 Quasar-Cluster Associations: Observations and Models 



clusters 


quasars 








{q)obs 


(ae/103)'^ 






ref 


Zwicky 


B < 18.5 


0.2 


1.8 


52 


1 7+0.5 
■■■•'-0.4 


r 0+1.6 
"3-'5_i.6 


28l?° 


o.ioiro^ 


1 


Abell 


5 > 2 Jy 


0.1 


2.0 


24 


1 7+0.5 
■■■•'-0.5 


4 7+1-2 

^•'-1.8 


99+13 
^^-14 


0.28t°:lo 


2 


UKJ287S 


B < 18.5 


0.15 


1.5 


7.2 


Z.U_Q 2 


9 q+0.2 
^■^-0.2 


^•"^-0.9 


12+0-02 

'-'•-'-^-0.02 


3 


Zwicky 


B <19 
5 > 1 Jy 


0.2 


1 


78 


~ 1.3 


4.3 
5.6 


18 
31 


0.06 
0.11 


4 



"Mean cluster redshift 
^Mcan quasar redshift 
'^Searcli range in arcminutes 

''Required cluster velocity dispersion in units of 1000 km/s 

*^Rcquircd surface mass density in g/cm^ for = 1 and Hq = 50 km/s/Mpc. 

-''References - (l)Rodrigues-Williams and Hogan, 1994; (2) Wu and Han, 1995; (3) Rodrigues and 
Hawkins, 1995; (4) Scitz and Schneider, 1995. 
9 Clusters in UKJ287 field 
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